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Abstract. Einstein's special theory of relativity revolutionized physics by 
teaching us that space and time are not separate entities, but join as "spacetime". 
His general theory of relativity further taught us that spacetime is not just a 
stage on which dynamics takes place, but is a participant: The field equation 
of general relativity connects matter dynamics to the curvature of spacetime. 
Curvature is responsible for gravity, carrying us beyond the Newtonian conception 
of gravity that had been in place for the previous two and a half centuries. Much 
research in gravitation since then has explored and clarified the consequences of 
this revolution; the notion of dynamical spacetime is now firmly established in 
the toolkit of modern physics. Indeed, this notion is so well established that we 
may now contemplate using spacetime as a tool for other science. One aspect of 
dynamical spacetime — its radiative character, "gravitational radiation" — will 
inaugurate entirely new techniques for observing violent astrophysical processes. 
Over the next one hundred years, much of this subject's excitement will come 
from learning how to exploit spacetime as a tool for astronomy. This article is 
intended as a tutorial in the basics of gravitational radiation physics. 



1. Introduction: Spacetime and gravitational waves 

Einstein's special relativity £Q taught us that space and time are not simply abstract, 
external concepts, but must in fact be considered measured observables, like any other 
quantity in physics. This reformulation enforced the philosophy that Newton sought 
to introduce in laying out his laws of mechanics [2]: 

... I frame no hypotheses; for whatever is not reduced from the phenomena is 
to be called an hypothesis; and hypotheses . . . have no place in experimental 
philosophy . . . 

Despite his intention to stick only with that which can be observed, Newton described 
space and time using exactly the abstract notions that he otherwise deplored 

Absolute space, in its own nature, without relation to anything external, 
remains always similar and immovable 

Absolute, true, and mathematical time, of itself, and from its own nature, 
flows equably without relation to anything external. 

Special relativity put an end to these abstractions: Time is nothing more than 
that which is measured by clocks, and space is that which is measured by rulers. The 
properties of space and time thus depend on the properties of clocks and rulers. The 
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constancy of the speed of light as measured by observers in different reference frames, 
as observed in the Michelson-Morley experiment, forces us inevitably to the fact that 
space and time are mixed into spacetime. Ten years after his paper on special relativity, 
Einstein endowed spacetime with curvature and made it dynamical 5 . This provided 
a covariant theory of gravity , in which all predictions for physical measurements are 
invariant under changes in coordinates. In this theory, general relativity, the notion 
of "gravitational force" is reinterpreted in terms of the behavior of geodesies in the 
curved manifold of spacetime. 

To be compatible with special relativity, gravity must be causal: Any change to 
a gravitating source must be communicated to distant observers no faster than the 
speed of light, c. This leads immediately to the idea that there must exist some notion 
of "gravitational radiation". As demonstrated by Bernard Schutz, one can actually 
calculate with surprising accuracy many of the properties of gravitational radiation 
simply by combining a time dependent Newtonian potential with special relativity [7] . 

The first calculation of gravitational radiation in general relativity is due to 
Einstein. His initial calculation jS] was "marred by an error in calculation" (Einstein's 
words), and was corrected in 1918 9^ (albeit with an overall factor of two error). 
Modulo a somewhat convoluted history (discussed in great detail by Kennefick |l()j ) 
owing (largely) to the difficulties of analyzing radiation in a nonlinear theory, Einstein's 
final result stands today as the leading-order "quadrupole formula" for gravitational 
wave emission. This formula plays a role in gravity theory analogous to the dipolc 
formula for electromagnetic radiation, showing that gravitational waves (hereafter 
abbreviated GWs) arise from accelerated masses exactly as electromagnetic waves 
arise from accelerated charges. 

The quadrupole formula tells us that GWs are difficult to produce — very large 
masses moving at relativistic speeds are needed. This follows from the weakness of 
the gravitational interaction. A consequence of this is that it is extremely unlikely 
there will ever be an interesting laboratory source of GWs. The only objects 
massive and relativistic enough to generate detectable GWs are astrophysical. Indeed, 
experimental confirmation of the existence of GWs has come from the study of binary 
neutron star systems — the variation of the mass quadrupole in such systems is 
large enough that GW emission changes the system's characteristics on a timescale 
short enough to be observed. The most celebrated example is the "Hulse- Taylor" 
pulsar, B1913+16, reported by Russell Hulse and Joe Taylor in 1975 ^5- Thirty 
years of observation have shown that the orbit is decaying; the results match with 
extraordinary precision general relativity's prediction for such a decay due to the 
loss of orbital energy and angular momentum by GWs. For a summary of the most 
recent data, see Ref. Fig- 1- Hulse and Taylor were awarded the Nobel Prize 
for this discovery in 1993. Since this pioneering system was discovered, several other 
double neutron star systems "tight" enough to exhibit strong GW emission have been 
discovered [T31 IT1 IT51 ITI)| . 

Studies of these systems prove beyond a reasonable doubt that GWs exist. What 
remains is to detect the waves directly and exploit them — to use GWs as a way 
to study astrophysical objects. The contribution to this volume by Danzmann 21 
discusses the challenges and program of directly measuring these waves. Intuitively, 
it is clear that measuring these waves must be difficult — the weakness of the 
gravitational interaction ensures that the response of any detector to gravitational 
waves is very small. Nonetheless, technology has brought us to the point where 
detectors are now beginning to set interesting upper limits on GWs from some sources 
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^lEHEBGII]- First direct detection is now hopefully not too far in the future. 

The real excitement will come when detection becomes routine. We will then 
have an opportunity to retool the "physics experiment" of direct detection into the 
development of astronomical observatories. Some of the articles appearing in this 
volume will discuss likely future revolutions which, at least conceptually, should change 
our notions of spacetime in a manner as fundamental as Einstein's works in 1905 and 
1915 (see, e.g., papers by Ashtekar and Horowitz). Such a revolution is unlikely to 
be driven by GW observations — most of what we expect to learn using GWs will 
apply to regions of spacetime that are well-described using classical general relativity; 
it is highly unlikely that Einstein's theory will need major revisions prompted by GW 
observations. Any revolution arising from GW science will instead be in astrophysics: 
Mature GW measurements have the potential to study regions of the Universe that 
are currently inaccessible to our instruments. During the next century of spacetime 
study, spacetime will itself be exploited to study our Universe. 

1.1. Why this article? 

As GW detectors have improved and approached maturity, many articles have been 
written reviewing this field and its promise. One might ask: Do we really need another 
one? As partial answer to this question, we note that Richard Price requested this 
article very nicely. More seriously, our goal is to provide a brief tutorial on the basics 
of GW science, rather than a comprehensive survey of the field. The reader who is 
interested in such a survey can find such in Refs. [^ll^3ll2105ll^Dll7fll^5109ll3l7| . 
Other reviews on the basics of GW science can be found in Refs. [21021; we also 
recommend the dedicated conference procedings [HHl Ell EH! • 

We assume that the reader has a basic familiarity with general relativity, at least 
at the level of Hartle's textbook [3^1; thus, we assume the reader understands metrics 
and is reasonably comfortable taking covariant derivatives. We adapt what Baumgarte 
and Shapiro |37| call the "Fortran convention" for indices: a, 6, c . . .h denote spacetime 
indices which run over 0, 1, 2, 3 or t, x, y, z, while i,j,k...n denote spatial indices which 
run over 1,2,3. We use the Einstein summation convention throughout — repeated 
adjacent indices in the "upstairs" and "downstairs" positions imply a sum: 



When we discuss linearized theory, we will sometimes be sloppy and sum over adjacent 
spatial indices in the same position. Hence, 



is valid in linearized theory. (As we will discuss in Sec. [3 this is allowable because, 
in linearized theory, the position of a spatial index is immaterial in Cartesian 
coordinates.) A quantity that is symmetrized on pairs of indices is written 



Throughout most of this article we use "relativist's units", in which G = 1 = c; mass, 
space, and time have the same units in this system. The following conversion factors 
are often useful for converting to "normal" units: 




a 




A(ab) = 7J (A a b + A ba ) . 



1 second = 299, 792, 458 meters ~ 3 x 10 s meters 
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1 Mq = 1476.63 meters ~ 1.5 kilometers 

= 4.92549 x 10~ 6 seconds ~ 5 microseconds . 

(1 Mq is one solar mass.) We occasionally restore factors of G and c to write certain 
formulae in normal units. 

Section |21 provides an introduction to linearized gravity, deriving the most basic 
properties of GWs. Our treatment in this section is mostly standard. One aspect of 
our treatment that is slightly unusual is that we introduce a gauge-invariant formalism 
that fully characterizes the linearized gravity's degrees of freedom. We demonstrate 
that the linearized Einstein equations can be written as 5 Poisson-type equations for 
certain combinations of the spacetime metric, plus a wave equation for the transverse- 
traceless components of the metric perturbation. This analysis helps to clarify which 
degrees of freedom in general relativity are radiative and which are not, a useful 
exercise for understanding spacetime dynamics. 

Section |31 analyses the interaction of GWs with detectors whose sizes are small 
compared to the wavelength of the GWs. This includes ground-based interferometric 
and resonant-mass detectors, but excludes space-based interferometric detectors. The 
analysis is carried out in two different gauges; identical results are obtained from both 
analyses. Section0|derives the leading-order formula for radiation from slowly-moving, 
weakly self-gravitating sources, the quadrupole formula discussed above. 

In Sec. we develop linearized theory on a curved background spacetime. Many 
of the results of "basic linearized theory" (Sec. |2J) carry over with slight modification. 
We introduce the "geometric optics" limit in this section, and sketch the derivation 
of the Isaacson stress-energy tensor, demonstrating how GWs carry energy and curve 
spacetime. Section provides a very brief synopsis of GW astronomy, leading the 
reader through a quick tour of the relevant frequency bands and anticipated sources. 
We conclude by discussing very briefly some topics that we could not cover in this 
article, with pointers to good reviews. 

2. The basic basics: Gravitational waves in linearized gravity 

The most natural starting point for any discussion of GWs is linearized gravity. 
Linearized gravity is an adequate approximation to general relativity when the 
spacetime metric, g a b, may be treated as deviating only slightly from a flat metric, 

Vab- 

gab = Vab + Kb, WhahW^l- (2.1) 

Here r\ a b is defined to be diag(— 1, 1,1,1) and \ \h a b\ \ means "the magnitude of a typical 
non-zero component of h a b" • Note that the condition ||/i a &|| "C 1 requires both the 
gravitational field to be weak, and in addition constrains the coordinate system to 
be approximately Cartesian. We will refer to h a b as the metric perturbation; as 
we will see, it encapsulates GWs, but contains additional, non-radiative degrees of 
freedom as well. In linearized gravity, the smallness of the perturbation means that 
we only keep terms which are linear in h a b — higher order terms are discarded. As 
a consequence, indices are raised and lowered using the flat metric T] a b- The metric 
perturbation h a b transforms as a tensor under Lorentz transformations, but not under 
general coordinate transformations. 
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We now compute all the quantities which are needed to describe linearized gravity. 
The components of the affine connection (Christoffel coefficients) are given by 

T a 6c = 7;V ad {d c hdb + dbhdc - ddhc) 

= \ (d c h a b + d b h\ - d a h bc ) . (2.2) 

Here d a means the partial derivative d/dx a . Since we use r\ a b to raise and lower 
indices, spatial indices can be written either in the "up" position or the "down" 
position without changing the value of a quantity: f x = f x . Raising or lowering a 
time index, by contrast, switches sign: / = —ft. The Riemann tensor we construct 
in linearized theory is then given by 

R a bcd — d c F a bd — ddF a bc 

= \ {d c d b h a d + d d d a h bc - dcd a h bd - d d d b h a c) . (2.3) 

From this, we construct the Ricci tensor 

Rab = R C acb = \ (d c d b h C a + d C d a h bc - Oh ab - d a d b h) , (2.4) 

where h = h a a is the trace of the metric perturbation, and □ = d c d c = V 2 — df is the 
wave operator. Contracting once more, we find the curvature scalar: 

R = R a a = (d c d a h c a - ah) (2.5) 

and finally build the Einstein tensor: 

G a b = Rab — -^HabR 

= 2 (d c dbh c a + d c d a h bc - Uh ab - d a d b h 

-Vabd c d d h c d + VabOh) . (2.6) 

This expression is a bit unwieldy. Somewhat remarkably, it can be cleaned up 
significantly by changing notation: Rather than working with the metric perturbation 
h ab , we use the trace-reversed perturbation h ab = h ab — \ij a bh. (Notice that h a a = —h, 
hence the name "trace reversed".) Replacing h ab with h ab + \r) a bh in Eq. (|2.6|l and 
expanding, we find that all terms with the trace h are canceled. What remains is 

Gab = \ (d c d b h c a + d c dahbc - ah ab - r ]ab d c d d h c <i ) . (2.7) 

This expression can be simplified further by choosing an appropriate coordinate 
system, or gauge. Gauge transformations in general relativity are just coordinate 
transformations. A general infinitesimal coordinate transformation can be written 
as x al = x a + £ a , where £ a (x b ) is an arbitrary infinitesimal vector field. This 
transformation changes the metric via 

h' ab - h ab - 29 (a 6) , (2.8) 

so that the trace-reversed metric becomes 

Kb = Kb - \^abti 

= h ab - 2d [b S, a) + T) ab d% . (2.9) 
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A class of gauges that are commonly used in studies of radiation are those satisfying 
the Lorentz gauge condition 

d a h ab = . (2.10) 

(Note the close analogy to Lorentz gaugef in electromagnetic theory, d a A a = 0, where 
A a is the potential vector.) 

Suppose that our metric perturbation is not in Lorentz gauge. What properties 
must £ a satisfy in order to impose Lorentz gauge? Our goal is to find a new metric 
h' ab such that d a h' ab = 0: 

d a h' ab = d a h ab ~ d a d b £ a - + d b 8% c (2.11) 
= 3 a h ab - . (2.12) 

Any metric perturbation h ab can therefore be put into a Lorentz gauge by making an 
infinitesimal coordinate transformation that satisfies 

□6 = d a h ab ■ (2.13) 

One can always find solutions to the wave equation i|2.13[l . thus achieving Lorentz 
gauge. The amount of gauge freedom has now been reduced from 4 freely specifiable 
functions of 4 variables to 4 functions of 4 variables that satisfy the homogeneous wave 
equation n£ 6 = 0, or, equivalently, to 8 freely specifiable functions of 3 variables on 
an initial data hypersurface. 

Applying the Lorentz gauge condition (|2.1U|) to the expression (|2.7|l for the 
Einstein tensor, we find that all but one term vanishes: 

G ab = —dhab. (2.14) 

Thus, in Lorentz gauges, the Einstein tensor simply reduces to the wave operator 
acting on the trace reversed metric perturbation (up to a factor —1/2). The linearized 
Einstein equation is therefore 

Uh ab = -16nT ab ; (2.15) 

in vacuum, this reduces to 

Uh ab = . (2.16) 

Just as in electromagnetism, the equation l|2.15|l admits a class of homogeneous 
solutions which are superpositions of plane waves: 

fc ot (x, t) = Re J d 3 k A ab (k)e^*-^ . (2.17) 

Here, lu = |k|. The complex coefficients A ab (k) depend on the wavevector k but are 
independent of x and t. They are subject to the constraint k a A ab = (which follows 
from the Lorentz gauge condition), with k a = (w,k), but are otherwise arbitrary. 
These solutions are gravitational waves. 

fFairly recently, it has become widely recognized that this gauge was in fact invented by Ludwig 
Lorenz, rather than by Hendrik Lorentz. The inclusion of the "t" seems most likely due to confusion 
between the similar names; see Ref. 1381 for detailed discussion. Following the practice of Griffiths 
f 1391 . p. 421), we bow to the weight of historical usage in order to avoid any possible confusion. 
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2.1. Globally vacuum spacetimes: Transverse traceless (TT) gauge 

We now specialize to globally vacuum spacetimes in which T a i, — everywhere, and 
which are asymptotically flat (for our purposes, h a b — > as r — > oo). Equivalently, we 
specialize to the space of homogeneous, asymptotically fiat solutions of the linearized 
Einstein equation (|2.15(l . For such spacetimes one can, along with choosing Lorentz 
gauge, further specialize the gauge to make the metric perturbation be purely spatial 

h tt =h u = (2.18) 

and traceless 

h = h^ = 0. (2.19) 

The Lorentz gauge condition l|2.10[l then implies that the spatial metric perturbation 
is transverse: 

dihij = 0. (2.20) 

This is called the transverse-traceless gauge, or TT gauge. A metric perturbation 
that has been put into TT gauge will be written ht^. Since it is traceless, there is no 
distinction between fij^ and h^. 

The conditions l|2.18|l and l|2.19|l comprise 5 constraints on the metric, while the 
residual gauge freedom in Lorentz gauge is parameterized by 4 functions that satisfy 
the wave equation. It is nevertheless possible to satisfy these conditions, essentially 
because the metric perturbation satisfies the linearized vacuum Einstein equation. 
When the TT gauge conditions are satisfied the gauge is completely fixed. 

One might wonder why we would choose TT gauge. It is certainly not necessary; 
however, it is extremely convenient, since the TT gauge conditions completely fix 
all the local gauge freedom. The metric perturbation h^ therefore contains only 
physical, non-gauge information about the radiation. In TT gauge there is a close 
relation between the metric perturbation and the linearized Riemann tensor R a bcd 
[which is invariant under the local gauge transformations (|2.8|l by Eq. 1)2. 3fl ] , namely 

1- 
-i 

2 

In a globally vacuum spacetime, all non-zero components of the Riemann tensor can 
be obtained from Ritjt via Ricmann's symmetries and the Bianchi identity. In a more 
general spacetime, there will be components that are not related to radiation; this 
point is discussed further in Sec. 12.21 

Transverse traceless gauge also exhibits the fact that gravitational waves have 
two polarization components. For example, consider a GW which propagates in the 
z direction: Kf^ = ft« T (f — z) is a valid solution to the wave equation C/iy T = 0. 
The Lorentz condition dji^j = implies that h^J{t — z) = constant. This constant 
must be zero to satisfy the condition that h a b — > as r — > oo. The only non-zero 
components of hj^ are then h^J, /iJJ, hyj, and hyj. Symmetry and the tracefree 
condition (|2.19(l further mandate that only two of these are independent: 

hJI= -h^ = h+(t-z); (2.22) 

h™ = h™ = h x {t-z). (2.23) 

The quantities h + and h x are the two independent waveforms of the GW (see Fig.^J. 

For globally vacuum spacetimes, one can always satisfy the TT gauge conditions. 
To see this, note that the most general gauge transformation £ a that preserves the 



Rim = -^W- ( 2 - 21 ) 
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Lorentz gauge condition (|2.1U|) satisfies D£ a = 0, from Eq. (|2.12|) . The general solution 
to this equation can be written 

£° = Re J d 3 k C Q (k)e i ( k x -" t ) (2.24) 

for some coefficients C°(k). Under this transformation the tensor A ab (k) in Eq. (|2.17Jl 
transforms as 

A ab -> A' ab = A ab - 2ik (a C b) + in ab k d C d . (2.25) 

Achieving the TT gauge conditions (|2.2U|) and (|2.19|l therefore requires finding, for 
each k, a C a (k) that satisfies the two equations 

= n ab A' ab = if h A ab + 2ik a C a (2.26) 
= A' tb = A tb - iC b kt - iC t k b + iS'hi^Ca) ; (2.27) 

is the Kronecker delta — zero for b ^ t, unity otherwise. An explicit solution to 
these equations is given by 

where k a = (w,k) and l a — (ui, — k). 



2.2. Global spacetimes with matter sources 

We now return to the more general and realistic situation in which the stress-energy 
tensor is non-zero. We continue to assume that the linearized Einstein equations are 
valid everywhere in spacetime, and that we consider asymptotically flat solutions only. 
In this context, the metric perturbation h ab contains (i) gauge degrees of freedom; 
(ii) physical, radiative degrees of freedom; and (iii) physical, non-radiative degrees 
of freedom tied to the matter sources. Because of the presence of the physical, 
non-radiative degrees of freedom, it is not possible in general to write the metric 
perturbation in TT gauge. However, the metric perturbation can be split up uniquely 
into various pieces that correspond to the degrees of freedom (i), (ii) and (iii), and 
the radiative degrees of freedom correspond to a piece of the metric perturbation that 
satisfies the TT gauge conditions, the so-called TT piece. 

This aspect of linearized theory is obscured by the standard, Lorentz gauge 
formulation (|2.15() of the linearized Einstein equations. There, all the components 
of h ab appear to be radiative, since all the components obey wave equations. In this 
subsection, we describe a formulation of linearized theory which focuses on gauge 
invariant observables. In particular, we will see that only the TT part of the metric 
obeys a wave equation in all gauges. We show that the non-TT parts of the metric 
can be gathered into a set of gauge invariant functions; these functions are governed 
by Poisson equations rather than wave equations. This shows that the non-TT pieces 
of the metric do not exhibit radiative degrees of freedom. Although one can always 
choose gauges like Lorentz gauge in which the non-radiative parts of the metric obey 
wave equations and thus appear to be radiative, this appearance is gauge artifact. 
Such gauge choices, although useful for calculations, can cause one to mistake purely 
gauge modes for truly physical radiation. 

Interestingly, the first analysis contrasting physical radiative degrees of freedom 
from purely coordinate modes appears to have been performed by Eddington in 1922 
40 . Eddington was somewhat suspicious of Einstein's analysis as Einstein chose 
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a gauge in which all metric functions propagated with the speed of light. Though the 
entire metric appeared to be radiative (by construction), Einstein found that only the 
"transverse-transverse" pieces of the metric carried energy. Eddington wrote: 

Weyl has classified plane gravitational waves into three types, viz. (1) 
longitudinal-longitudinal; (2) longitudinal-transverse; (3) transverse-trans- 
verse. The present investigation leads to the conclusion that transverse- 
transverse waves are propagated with the speed of light in all systems of 
co-ordinates. Waves of the first and second types have no fixed velocity 
- a result which rouses suspicion as to their objective existence. Einstein 
had also become suspicious of these waves (in so far as they occur in his 
special co-ordinate system) for another reason, because he found that they 
convey no energy. They are not objective, and (like absolute velocity) are 
not detectable by any conceivable experiment. They are merely sinuosities in 
the co-ordinate system, and the only speed of propagation relevant to them 
is "the speed of thought." 

... It is evidently a great convenience in analysis to have all waves, both 
physical and spurious, travelling with one velocity; but it is liable to obscure 
physical ideas by mixing them up so completely. The chief new point in 
the present discussion is that when unrestricted co-ordinates are allowed the 
genuine waves continue to travel with the velocity of light and the spurious 
waves cease to have any fixed velocity. 

Unfortunately, Eddington's wry dismissal of unphysical modes as propagating with 
"the speed of thought" is often taken by skeptics (and crackpots) as applying to all 
gravitational perturbations. Eddington in fact showed quite the opposite. We do so 
now using somewhat more modern notation; our presentation is essentially the flat- 
spacetime limit of Bardeen's |41| gauge-invariant cosmological perturbation formalism. 
A similar treatment can be found in lecture notes by Bertschinger |42| . 

We begin by defining the decomposition of the metric perturbation h a b, in any 
gauge, into a number of irreducible pieces. Assuming that h a b — ► as r — * oo, we 
define the quantities 0, 7, H, Si, A and hj^ via the equations 

h tt = 20 , (2.29) 
hu = fir + 5 l7 , (2.30) 

hj = hf? + ^HSij + + (didj - ISijV 2 ) A , (2.31) 

together with the constraints 

difSi =0 (1 constraint) (2.32) 

diEi =0 (1 constraint) (2.33) 

dihjf = (3 constraints) (2.34) 

5 ij h^ = (1 constraint) (2.35) 

and boundary conditions 

7^0, Si -> 0, A -► 0, V 2 A -> (2.36) 

as r — > 00. Here H = S lJ hij is the trace of the spatial portion of the metric 
perturbation, not to be confused with the spacetime trace h = r] ab h a b that we used 
earlier. The spatial tensor hj^- is transverse and traceless, and is the TT piece of 
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the metric discussed above which contains the physical radiative degrees of freedom. 
The quantities fa and 9,7 are the transverse and longitudinal pieces of h t i- The 
uniqueness of this decomposition follows from taking a divergence of Eq. (|2.30|) giving 
V 2 7 = dih t i, which has a unique solution by the boundary condition (|2.36(1 . Similarly 
taking two derivatives of Eq. I|2.31[l yields the equation 2V 2 V 2 A = 3didjhij — V 2 iJ, 
which has a unique solution by Eq. (|2.36|1 . Having solved for A, one can obtain a 
unique £$ by solving 3V 2 £i = 6djhij — 2diH — idi V 2 A. 

The total number of free functions in the parameterization 12.29JI - (|2.31|1 of the 
metric is 16: 4 scalars (</>, 7, H, and A), 6 vector components (fa and and 6 
symmetric tensor components (hjj). The number of constraints (|2.32|l - l|2.35|l is 6, 
so the number of independent variables in the parameterization is 10, consistent with 
a symmetric 4x4 tensor. 

We next discuss how the variables <f>, fa, 7, H, Ei, A and hJ } T transform under 
gauge transformations £ a with £ a — > as r — > 00. We parameterize such gauge 
transformation as 

& = = (A, Bi + (2.37) 

where diBi = and C — > as r — > 00; thus Bi and diC are the transverse and 
longitudinal pieces of the spatial gauge transformation. The transformed metric is 
h a b — 2<9( £ f) ); decomposing this transformed metric into its irreducible pieces yields 
the transformation laws 

(2.38) 
(2.39) 
(2.40) 
(2.41) 
(2.42) 
(2.43) 
(2.44) 

Gathering terms, we see that the following combinations of these functions are gauge 
invariant: 

$ ee -0 + 7-Ia, (2.45) 
9 ee i (H - V 2 A) , (2.46) 

= fa - ^ ; (2.47) 

hj? is gauge-invariant without any further manipulation. In the Newtonian limit $ 
reduces to the Newtonian potential $_/v, while 6 = — 2$jv- The total number of free, 
gauge-invariant functions is 6: 1 function 0; 1 function $; 3 functions 3$, minus 1 
due to the constraint <9jSj = 0; and 6 functions hj? , minus 3 due to the constraints 
dihj^- = 0, minus 1 due to the constraint 8 V hj? = 0. This is in keeping with the 
fact that in general the 10 metric functions contain 6 physical and 4 gauge degrees of 
freedom. 

We would now like to enforce Einstein's equation. Before doing so, it is useful 
to first decompose the stress energy tensor in a manner similar to that of our 





4>-A, 


fa - 


fa-Bi, 


7 -t 


1-A-C, 


H -> 


H-2V 2 C , 


A -> 


A-2C , 


e« -» 


£< — 25j , 




n v ■ 
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decomposition of the metric. We define the quantities p, Si, S, P, dy, cr, and a 
via the equations 

T tt =p, (2.48) 

T u =Si + dtS , (2.49) 



T ti = PSij + o-ij + d^o-j) + didj - V 2 ) a, (2.50) 



together with the constraints 





= 0, 


(2.51) 


BiCFi 


= 0, 


(2.52) 


di(Tij 


= 0, 


(2.53) 




= 0, 


(2.54) 



and boundary conditions 

S^0, cr,- -*•(), cr^O, V 2 cr^0 (2.55) 

as r — > co. These quantities are not all independent. The variables p, P, Si and o~ij 

can be specified arbitrarily; stress-energy conservation (d a T a i, = 0) then determines 
the remaining variables S, a, and <Xj via 

V 2 S = p , (2.56) 

V 2 a = ~\P+\S, (2.57) 

V 2 fi 4 = 2S t . (2.58) 

We now compute the Einstein tensor from the metric (|2.29|l - H2.31|) . The result 
can be expressed in terms of the gauge invariant observables: 

G tt = - V 2 9 , (2.59) 

G tl = - Jv 2 S 4 - d t O , (2.60) 



2 

-v 2 (2$ + e)-e 



Gn = - ^OhJ? - d(iEj) - ^didj (2$ + 9) 



2 

(2.61) 

We finally enforce Einstein's equation G a b — SnTa^ and simplify using the conservation 
relations H2.56[) ~ (|2.58l) : this leads to the following field equations: 

V 2 9 = -87175, (2.62) 
V 2 $ = 4?r (p + 3P - 35) , (2.63) 



V 2 S, = - 167r5i , (2.64) 
z^ T = - IQno-ij 



□/i,TT = _ 16 _ (2.65) 



Notice that only the metric components hj^ obey a wave-like equation. 

The other variables 0, and 3j are determined by Poisson-type equations. Indeed, in 
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a purely vacuum spacetime, the field equations reduce to five Laplace equations and 
a wave equation: 



v 2 e vac 


= o, 


(2.66) 


v 2 $ vac 


= o, 


(2.67) 


vac 


= o, 


(2.68) 


hl TT.vac 


= o . 


(2.69) 



This manifestly demonstrates that only the metric components — the transverse, 
traceless degrees of freedom of the metric perturbation — characterize the radiative 
degrees of freedom in the spacetime. Although it is possible to pick a gauge in which 
other metric components appear to be radiative, they will not be: Their radiative 
character is an illusion arising due to the choice of gauge or coordinates. 

The field equations (|2.62|) - 1)2. 65|) also demonstrate that, far from a dynamic, 
radiating source, the time-varying portion of the physical degrees of freedom in the 
metric is dominated by hj^ . If we expand the gauge invariant fields $, 6, Si and 
hj? in powers of 1/r, then, at sufficiently large distances, the leading-order 0(l/r) 
terms will dominate. For the fields 0, $ and Hj, the coefficients of the 1/r pieces are 
simply the conserved mass J d 3 xp or the conserved linear momentum — J cPxSi, from 
the conservation relations (|2.56|) - (|2.58|) . Thus, the only time- varying piece of the 
physical degrees of freedom in the metric perturbation at order 0(1/ r) is the TT piece 
hj^ . An alternative proof of this result is given in Exercise 19.1 of Misner, Thorne 
and Wheeler 0]. 

Although the variables $, 8, and hJ 3 T have the advantage of being gauge 
invariant, they have the disadvantage of being non-local. Computation of these 
variables at a point requires knowledge of the metric perturbation h a {> everywhere. 
This non-locality obscures the fact that the physical, non-radiative degrees of freedom 
are causal, a fact which is explicit in Lorentz gauge f. On the other hand, many 
observations that seek to detect GWs are sensitive only to the value of the Riemann 
tensor at a given point in space (see Sec. |3J). For example, the Riemann tensor 
components Ritjt, which are directly observable by detectors such as LIGO, are given 
in terms of the gauge invariant variables as 

Ritjt = -\hj? + + S (w) - ^6%. (2.70) 

Thus, at least certain combinations of the gauge invariant variables are locally 
observable. 

2. 3. Local regions of spacetime 

In the previous subsection we described a splitting of metric perturbations into 
radiative, non-radiative, and gauge pieces. This splitting requires that the linearized 
Einstein equations be valid throughout the spacetime. However, this assumption is 
not valid in the real Universe: Many sources of GWs are intrinsically strong field 
sources and cannot be described using linearized theory, and on cosmological scales 
the metric of our Universe is not close to the Minkowski metric. Furthermore the 
splitting requires a knowledge of the metric throughout all of spacetime, whereas 

fOne way to see that the guage invariant degrees of freedom are causal is to combine the vacuum 
wave equation 12.161 for the metric perturbation with the expression 12.31 for the gauge-invariant 
Riemann tensor. This gives the wave equation Dfl^j = 0. 
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any measurements or observations can probe only finite regions of spacetime. For 
these reasons it is useful to consider linearized perturbation theory in finite regions of 
spacetime, and to try to define gravitational radiation in this more general context. 

Consider therefore a finite volume V in space. Can we split up the metric 
perturbation h a b in V into radiative and non-radiative pieces? In general, the answer 
is no: Within any finite region GWs cannot be distinguished from time- varying near- 
zone fields generated by sources outside that region. One way to see this is to note that 
in finite regions of space, the decomposition of the metric into various pieces becomes 
non-unique, as does the decomposition of vectors into transverse and longitudinal 
pieces. [For example the vector (a; 2 — y 2 )d z is both transverse and longitudinal.] 
Alternatively, we note that within any finite vacuum region V, one can always find a 
gauge which is locally TT, that is, a gauge which satisfies the conditions (|2.18J) - (|2.2()|l 
within the region. (This fact does not seem to be widely known, so we give a proof in 
|Appendix A| ). In particular, this applies to the static Coulomb-type field of a point 
source, as long as the source itself is outside of V. Consequently, isolating the TT piece 
of the metric perturbation does not yield just the radiative degrees of freedom within a 
local region - a TT metric perturbation may also contain, for example, Coulomb-type 
fields. 

Within finite regions of space, therefore, GWs cannot be defined at a fundamental 
level - one simply has time-varying gravitational fields. However, there is a certain 
limit in which GWs can be approximately defined in local regions, namely the limit 
in which the wavelength of the waves is much smaller than length and time scales 
characterizing the background metric. This definition of gravitational radiation is 
discussed in detail and in a more general context in Sec. 03 As discussed in that 
section, this limit will always be valid when one is sufficiently far from all radiating 
sources. 



3. Interaction of gravitational waves with a detector 

The usual notion of "gravitational force" disappears in general relativity, replaced 
instead by the idea that freely falling bodies follow geodesies in spacetime. Given a 
spacetime metric g a b and a set of spacetime coordinates x a , geodesic trajectories are 
given by the equation 

d 2 x a „ dx b dx c , . 

where r is proper time as measured by an observer travelling along the geodesic. By 
writing the derivatives in the geodesic equation (|3.1|l in terms of coordinate time t 
rather than proper time r, and by combining the a = t equation with the spatial, 
a = j equations, we obtain an equation for the coordinate acceleration: 

-(It + 2IV + + u*(I* + 2IV + T > 3 v k ), (3.2) 



—jjj- — tf-r ^ tj v ' jk ) ' \ tt ' ij x jk l 

where v % = dx l jdt is the coordinate velocity. 

Let us now specialize to linearized theory, with the non-flat part of our metric 
dominated by a GW in TT gauge. Further, let us specialize to non-relativistic motion 
for our test body. This implies that v l -C 1, and to a good approximation we can 
neglect the velocity dependent terms in Eq. I|3.2|l : 

d 2 x i 

-ar+i*— 0. (3-3) 
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In linearized theory and TT gauge, 

r*« = r m = \ (2d t h] t T - djhY) = o (3.4) 

since ft,J t T = 0. Hence, we find that d 2 x l /dt 2 = 0. 

Does this result mean that the GW has no effect? Certainly not! It just tells 
us that in TT gauge the coordinate location of a slowly moving, freely falling body is 
unaffected by the GW. In essence, the coordinates move with the waves. 

This result illustrates why, in general relativity, it is important to focus upon 
coordinate-invariant observables — a naive interpretation of the above result would 
be that freely falling bodies are not influenced by GWs. In fact the GWs cause the 
proper separation between two freely falling particles to oscillate, even if the coordinate 
separation is constant. Consider two spatial freely falling particles, located at z = 0, 
and separated on the x axis by a coordinate distance L c . Consider a GW in TT gauge 
that propagates down the z axis, hj^-(t, z). The proper distance L between the two 
particles in the presence of the GW is given by 



L 



dx y/g^ = I dx \Jl + ti££(t, z = 0) 



o Jo 

(3.5) 



dx 



1 + ^L T (M = 0) 



= L, 



1 + ^J(M = 0) 



Notice that we use the fact that the coordinate location of each particle is fixed in 
TT gauge! In a gauge in which the particles move with respect to the coordinates, 
the limits of integration would have to vary. Equation <|3.5|) tells us that the proper 
separation of the two particles oscillates with a fractional length change SL/L given 

by 

f 4^(M = 0). (3.6) 

Although we used TT gauge to perform this calculation, the result is gauge 
independent; we will derive it in a different gauge momentarily. Notice that h^J 
acts as a strain — a fractional length change. The magnitude h of a wave is often 
referred to as the "wave strain" . The proper distance we have calculated here is 
a particularly important quantity since it directly relates to the accumulated phase 
which is measured by laser interferometric GW observatories (cf. the contribution by 
Danzmann in this volume). The "extra" phase 5<j) accumulated by a photon that 
travels down and back the arm of a laser interferometer in the presence of a GW is 
S(f) — AttSL/X, where A is the photon's wavelength and 5L is the distance the end 
mirror moves relative to the beam splitterf . We now give a different derivation of the 
fractional length change l|3.6[) based on the concept of geodesic deviation. Consider 
a geodesic in spacetime given by x a — z a (r), where r is the proper time, with four 
velocity u°(r) = dz a /dr. Suppose we have a nearby geodesic x a (r) — z a {r) + L a (r), 
where L a (r) is small. We can regard the coordinate displacement L a as a vector 
L = L a d a on the first geodesic; this is valid to first order in L. Without loss of 
generality, we can make the connecting vector be purely spatial: L a u a = 0. Spacetime 
curvature causes the separation vector to change with time — the geodesies will move 

fThis description of the phase shift only holds if L <g A, so that the metric perturbation does 
not change value very much during a light travel time. This condition will be violated in the high 
frequency regime for space-based GW detectors; a more careful analysis of the phase shift is needed 
in this case 1431 . 
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further apart or closer together, with an acceleration given by the geodesic deviation 
equation 

u b V b {u c V c L a ) = -R\ cd [z(r)]u b L c u d ; (3.7) 

see, e.g., Ref. Chap. 21. This equation is valid to linear order in L a ; fractional 
corrections to this equation will scale as L/C, where £ is the lengthscale over which 
the curvature varies. 

For application to GW detectors, the shortest such lengthscale C is the wavelength 
A of the GWs. Thus, the geodesic deviation equation will have fractional corrections 
of order L/X. For ground-based detectors L <, & few km, while A <; 3000km (see Sec. 
16.111 : thus the approximation will be valid. For detectors with L <; A (e.g. the space 
based detector LISA) the analysis here is not valid and other techniques must be used 
to analyze the detector. 

A convenient coordinate system for analyzing the geodesic deviation equation 
l|3.7|) is the local proper reference frame of the observer who travels along the first 
geodesic. This coordinate system is defined by the requirements 

2 l (r) = 0, g ab (t,0) = Vab, r£.(t,0)=0, (3.8) 

which imply that the metric has the form 

ds 2 = -dt 2 + rfx 2 + O (^] . (3.9) 



Here 1Z is the radius of curvature of spacetime, given by 1Z~ 2 ~ ||-R a i>cd||- It also 
follows from the gauge conditions l|3.8|l that proper time r and coordinate time t 
coincide along the first geodesic, that u = d t and that L a = (0, L l ). 

Consider now the proper distance between the two geodesies, which are located 
at x 1 — and x % = L' l (t). From the metric (|3.9|l we see that this proper distance is just 
|L| = y/LiLi, up to fractional corrections of order L 2 /1Z 2 . For a GW of amplitude h 
and wavelength A we have ~Rr 2 ~ h/X 2 , so the fractional errors are ~ hL 2 /X 2 . (Notice 
that 1Z r-j C/Vh — the wave's curvature scale 1Z is much larger than the lengthscale 
C characterizing its variations.) Since we are restricting attention to detectors with 
L <^ A, these fractional errors are much smaller than the fractional distance changes 
~ h caused by the GW [Eq. (|3.6J) ]. Therefore, we can simply identify |L| as the proper 
separation. 

We now evaluate the geodesic deviation equation (|3.7J) in the local proper reference 
frame coordinates. From the conditions l|3.8|) it follows that we can replace the 
covariant time derivative operator u a V a with d/(dt). Using u = dt and L a = (0,L l ) 
we get 

d 2 LHt) , , • , , 

-^l = -R itjt (t,0)L : >(t). (3.10) 

Note that the key quantity entering into the equation, Rnjt, is gauge invariant in 
linearized theory, so we can use any convenient coordinate system to evaluate it. 
Using the expression (|2.21(l for the Riemann tensor in terms of the TT gauge metric 
perturbation hj^ we find 

dt* "2 dt 2 (6 J 

Integrating this equation using L l (f) = L % + 8L l (t) with |<5L| <C |Lo| gives 

5U{t) = \hJ?{t)L{. (3.12) 
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Figure 1. Lines of force for a purely + GW (left), and for a purely X GW (right). 
Figure kindly provided by Kip Thorne; originally published in Ref. 1441 . 

This equation is ideal for analyzing an interferometric GW detector. We choose 
Cartesian coordinates such that the interferometer's two arms lie along the x and y 
axes, with the beam splitter at the origin. For concreteness, let us imagine that the 
GW propagates down the z-axis. Then, as discussed in Sec. 12.11 the only non-zero 
components of the metric perturbation are ti^J = —hyy = h+ and hjy = hyj — h x , 
where h+ (t — z) and h x (t — z) are the two polarization components. We take the ends 
of one of the interferometer's two arms as defining the two nearby geodesies; the first 
geodesic is defined by the beam splitter at x = 0, the second by the end-mirror. From 
Eq. i|3.12f> we then find that the distances L = |L| of the arms' ends from the beam 
splitter vary with time as 



SL X 




L 


2 + ' 


SLy 






= -h 


L 


2 



(3.13) 



(Here the subscripts x and y denote the two different arms, not the components 
of a vector). These distance changes are then measured via laser interferometry. 
Notice that the GW (which is typically a sinusoidally varying function) acts tidally, 
squeezing along one axis and stretching along the other. In this configuration the 
detector is sensitive only to the + polarization of the GW. The x polarization acts 
similarly, except that it squeezes and stretches along a set of axes that are rotated 
with respect to the x and y axes by 45°. The force lines corresponding to the two 
different polarizations are illustrated in Fig. ^ 

Of course, we don't expect nature to provide GWs that so perfectly align 
with our detectors. In general, we will need to account for the detector's antenna 
pattern, meaning that we will be sensitive to some weighted combination of the two 
polarizations, with the weights depending upon the location of a source on the sky, 
and the relative orientation of the source and the detector. See Ref. j32J, Eqs. (104a, b) 
and associated text for further discussion. 

Finally, in our analysis so far of detection we have assumed that the only 
contribution to the metric perturbation is the GW contribution. However, in reality 
time- varying near zone gravitational fields produced by sources in the vicinity of the 
detector will also be present. From Eq. H3.10fl we see that the quantity that is actually 
measured by interferometric detectors is the space-time-space-time or electric-type 
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piece Rujt of the Riemann tensor (or more precisely the time-varying piece of this 
within the frequency band of the detector). From the general expression (|2.7U|) for 
this quantity we see that Rujt contains contributions from both hj^ describing GWs, 
and also additional terms describing the time-varying near zone gravitational fields. 
There is no way for the detector to separate these two contributions, and the time- 
varying near zone gravitational fields produced by motions of bedrock, air, human 
bodies, and tumbleweeds can all contribute to the output of the detector and act as 
sources of noise |47fl 14*71 HB] . 



4. The generation of gravitational waves: Putting in the source 

4-1- Slow-motion sources in linearized gravity 

Gravitational waves are generated by the matter source term on the right hand side 
of the linearized Einstein equation 

Uhab = -16^T afc , (4.1) 

cf. Eq. I|2.15|> (presented here in Lorentz gauge). In this section we will compute the 
leading order contribution to the spatial components of the metric perturbation for a 
source whose internal motions are slow compared to the speed of light ( "slow-motion 
sources"). We will then compute the TT piece of the metric perturbation to obtain 
the standard quadrupolc formula for the emitted radiation. 

Equation (|4.1I) can be solved by using a Green's function. A wave equation with 
source generically takes the form 

□/(t,x)=*(t,x), (4.2) 

where f(t, x) is the radiative field, depending on time t and position x, and s(t, x) is 
a source function. The Green's function G(t, x; t' , x') is the field which arises due to 
a delta function source; it tells how much field is generated at the "field point" (t, x) 
per unit source at the "source point" (i',x'): 

DG{t, x; t', x') = 5(t - t')6{x - x') . (4.3) 

The field which arises from our actual source is then given by integrating the Green's 
function against s(t, x): 

/(t,x) = J dt'd 3 x' G(t,x;t',x')s(t',x') . (4.4) 

The Green's function associated with the wave operator □ is very well known (see, 
e.g., B»: 

ai^w— l*-*-:*-?™ . (45) 

47r|x — x'| 

The quantity t— |x — x'|/c is the retarded time; it takes into account the lag associated 
with the propagation of information from events at x to position x'. The speed of 
light c has been restored here to emphasize the causal nature of this Green's function; 
we set it back to unity in what follows. 

Applying this result to Eq. we find 

m.,,)^/^ /-"; 1 -; 1 '" (4.6) 
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As already mentioned, the radiative degrees of freedom are contained entirely in the 
spatial part of the metric, projected transverse and traceless. First, consider the 
spatial part of the metric: 



, 3 ,T^-|x-x'|,x') 



Mi,x) = 4 / dV i — =— ! -i-U— i . (4.7) 



We have raised indices on the right-hand side, using the rule that the position of 
spatial indices in linearized theory is irrelevant. 

We now evaluate this quantity at large distances from the source. This allows 
us to replace the factor |x — x'| in the denominator with r = |x|. The corresponding 
fractional errors scale as ~ L/r, where L is the size of the source; these errors can be 
neglected. We also make the same replacement in the time argument of T,j: 

T l3 (t - |x - x'|,x') a T«(i - r, x'). (4.8) 

Using the formula |x — x'| = r — n l x 1 + 0(l/r) where n l = x t /r, we see that the 
fractional errors generated by the replacement l|4. 8J) scale as L/t, where r is the 
timescale over which the system is changing. This quantity is just the velocity of 
internal motions of the source (in units with c = 1), and is therefore small compared 
to one by our assumption. These replacements give 

hij (t, x) = - J d 3 x' T lj (t - r, x') , (4.9) 

which is the first term in a multipolar expansion of the radiation field. 

Equation (|4.9() almost gives us the quadrupole formula that describes GW 
emission (at leading order). To get the rest of the way there, we need to massage this 
equation a bit. The stress-energy tensor must be conserved, which means d a T ah = 
in linearized theory. Breaking this up into time and space components, we have 

d t T tt + d l T tl =0, (4.10) 

d t T u + djT ij = . (4.11) 
From this, it follows rather simply that 

dfT tt = d k diT kl . (4.12) 
Multiply both sides of this equation by x % x 3 . We first manipulate the left-hand side: 

d 2 t T u x l x 3 = df (T^V) . (4.13) 
Next, manipulate the right-hand side of Eq. (|4.12l) . multiplied by x % x 3 : 

d k d l T kl x i x j = d k di (T kl x l x J ) - 2d k (T ik x J + T kj x l ) + 2T lj . (4.14) 

This identity is easily verifiedf by expanding the derivatives and applying the identity 
diX 3 = Si 3 . We thus have 

dl (T u x l x j ) = d k di (T w a:V) - 2d k (T ik x j + T kj x l ) + 2T l] . (4.15) 

fAlthough one of us (SAH) was unable to do this simple calculation while delivering lectures at a 
summer school in Brownsville, TX. Never attempt to derive the quadrupole formula while medicated. 
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This yields 



d 3 x' T;. 



d 3 x' 



-df (T tt x H x' j ) + d k (T tk x' j + T kj x H ) 



-d k dt (T kl x H x'i) 



d 3 x'd 2 t (T tt x H x' j ) 



r ot z 

2 -^ld*x' P x'W 



(4.16) 



In going from the first line to the second, we used the fact that the second and third 
terms under the integral are divergences. Using Gauss's theorem, they can thus be 
recast as surface integrals; taking the surface outside the source, their contribution is 
zero. In going from the second line to the third, we used the fact that the integration 
domain is not time dependent, so we can take the derivatives out of the integral. 
Finally, we used the fact that T u is the mass density p. Defining the second moment 



Iij of the mass distribution via 



I l0 (t) = J d 3 x' p{t,x!)x n x' J , 

and combining Eqs. I|4.9(l and l|4.16|l now gives 

- 2 dH ij {t-r) 
lAh j "r dt 2 

When we subtract the trace from Zy, we obtain the quadrupole moment tensor: 



(4.17) 



(4.18) 



(4.19) 



This tensor will prove handy shortly. 

To complete the derivation, we must project out the non-TT pieces of the right- 
hand side of Eq. I|4.18|) . Since we are working to leading order in 1/r, at each field 
point x this operation reduces to algebraically projecting the tensor perpendicularly 
to the local direction of propagation n = x/r, and subtracting off the trace. It is 
useful to introduce the projection tensor, 

Pij = % - riinj . (4.20) 

This tensor eliminates vector components parallel to n, leaving only transverse 
components. Thus, 

hT 



"ij — hklPikPjl 

is a transverse tensor. Finally, we remove the trace; what remains is 
hklPikPjl - 7;PijPkihki ■ 



h TT 



Substituting Eq. (|4~T%t into (|4~221) . we obtain our final quadrupole formula: 



x) 



2d 2 I kl (t- 
r dt 2 



P ik {n)P jl {n)--P u (n)P i M) 



(4.21) 
(4.22) 

(4.23) 
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4-2. Extension to sources with non-negligible self gravity 

Our derivation of the quadrupole formula l|4.23|) assumed the validity of the linearized 
Einstein equations. In particular, the derivation is not applicable to systems with 
weak (Newtonian) gravity whose dynamics are dominated by self-gravity, such as 
binary star systemsf. This shortcoming of the above lincarizcd-gravity derivation 
of the quadrupole formula was first pointed out by Eddington. However, it is very 
straightforward to extend the derivation to encompass systems with non-negligible self 
gravity. 

In full general relativity we define the quantity h ab via 

^J~gg ab = Tf h - h ab , (4.24) 

where rf h ee diag (— 1, 1, 1, 1). When gravity is weak this definition coincides with our 
previous definition of h ab as a trace-reversed metric perturbation. We impose the 
harmonic gauge condition 

da(V~99 ab ) = d a h ab = 0. (4.25) 

In this gauge the Einstein equation can be written 

U &at h ab = -167r(T ab + t ab ) (4.26) 

where n aat = ri ab d a db is the flat-spacetime wave operator, and t ab is a pseudotensor 
that is constructed from h ab . Taking a coordinate divergence of this equation and 
using the gauge condition lj4.25|) shows that stress-energy conservation can be written 

d a {T ab + t ab ) = 0. (4.27) 

Equations H4.25(l . (|4.26() and 14.2711 are precisely the same equations as are used 
in the linearized-gravity derivation of the quadrupole formula, except for the fact that 
the stress energy tensor T ab is replaced by T ab + t ab . Therefore the derivation of the 
last subsection carries over, with the one modification that the formula (|4.17|) for ijj 
is replaced by 

I l0 (t) = J d 3 x' [T"(t,x') + t"(t,x')] x' l x' 3 . (4.28) 

In this equation the term t describes gravitational binding energy, roughly speaking. 
For systems with weak gravity, this term is negligible in comparison with the term 
T u describing the rest-masses of the bodies. Therefore the quadrupole formula (|4.23|l 
and the original definition (|4. 1 T|l of Itj continue to apply to the more general situation 
considered here. 



4-3. Dimensional analysis 

The rough form of the leading GW field that we just derived, Eq. (|4.23(l . can be 
deduced using simple physical arguments. First, we define some moments of the mass 
distribution. The zcroth moment is just the mass itself: 

M ee J pd 3 x = M . (4.29) 

fStress energy conservation in linearized gravity, d a T a i, = 0, forces all bodies to move on geodesies 
of the Minkowski metric. 
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(More accurately, this is the total mass-energy of the source.) Next, we define the 
dipole moment: 

Mi = J p Xl d 3 x = ML, . (4.30) 

Li is a vector with the dimension of length; it describes the displacement of the center 
of mass from our chosen origin. (As such, Mi is clearly not a very meaningful quantity 
— we can change its value simply by choosing a different origin.) 

If our mass distribution exhibits internal motion, then moments of the moss 
current, ji = pvi, are also important. The first moment is the spin angular momentum: 



Si = J pvj x k tijk d x — Si . (4-31) 

Finally, we look at the second moment of the mass distribution: 

M 2 = J px, Xj d 3 x = MLij (4.32) 

where Lij is a tensor with the dimension length squared. 

Using dimensional analysis and simple physical arguments, it is simple to see that 
the first moment that can contribute to GW emission is Consider first Mq. We 
want to combine Mq with the distance to our source, r, in such a way as to produce a 
dimensionless wavestrain h. The only way to do this (bearing in mind that the strain 
should fall off as 1/r, and restoring factors of G and c) is to put 

h~%^. (4.33) 

cr r 

Does this formula make sense for radiation? Not at all! Conservation of mass-energy 
tells us that Mq for an isolated source cannot vary dynamically. This h cannot be 
radiative; it corresponds to a Newtonian potential, rather than a GW. 

How about the moment Mi? In order to get the dimensions right, we must take 
one time derivative: 

(The extra factor of c converts the dimension of the time derivative to space, so that 
the whole expression is dimensionless.) Think carefully about the derivative of Ml: 

*¥± = ±Jp Xi <Px=[pv i <Px = P i . (4.35) 

This is the total momentum of our source. Our guess for the form of a wave 
corresponding to Mi becomes 
C P 

h~±-. (4.36) 

cr r 

Can this describe a GW? Again, not a chance: The momentum of an isolated source 
must be conserved. By boosting into a different Lorentz frame, we can always set 
P = 0. Terms like this can only be gauge artifacts; they do not correspond to 
radiation. [Indeed, terms like (|4.36l) appear in the metric of a moving black hole, 
and correspond to the relative velocity of the hole and the observer; see [5U], Chapter 
5-] 

How about Si? Dimensional analysis tells us that radiation from Si must take 
the form 

G d S\ . 
h~—- — -. (4.37) 

c 4 dt r y ' 
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Conservation of angular momentum tells us that the total spin of an isolated system 
cannot change, so we reject this term for the same reason that we rejected l|4.33[l - 
it cannot correspond to radiation. 
Finally, we examine M%: 

GcP_M2 

There is no conservation principle that allows us to reject this term. Comparing to Eq. 
H4.23f) . we see that this is the quadrupole formula we derived earlier, up to numerical 
factors. 

In "normal" units, the prefactor of this formula turns out to be G/c 4 — a small 
number divided by a very big number. In order to generate interesting amounts of 
GWs, the quadrupole moment's variation must be enormous. The only interesting 
sources of GWs will be those which have very large masses undergoing extremely 
rapid variation; even in this case the strain we expect from typical sources is tiny. 
The smallness of GWs reflects the fact that gravity is the weakest of the fundamental 
interactions. 



4-4 ■ Numerical estimates 

Consider a binary star system, with stars of mass mi and mi in a circular orbit with 
separation R. The quadrupole moment is given by 



■^ij — ( X%Xj 



(4.39) 



where [i = m\m2/(mi + m-i) is the binary's reduced mass and x is the relative 
displacement, with |x| = R. We use the center-of-mass reference frame, and choose 
the coordinate axes so that the binary lies in the xy plane, so x — x\ = Rcosttt, 
y = X2 = R sin fit, z — x% — 0. Let us further choose to evaluate the field on the z 
axis, so that n points in the z-direction. The projection operators in Eq. (|4.23|l then 
simply serve to remove the zj components of the tensor. Bearing this in mind, the 
quadrupole formula (|4.23l) yields 



■ TT 
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The quadrupole moment tensor is 
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its second derivative is 
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The magnitude h of a typical non-zero component of hlj is 
, 4[in 2 R 2 4 A iAf 2 / 3 ^ 2 / 3 



(4.40) 



(4.41) 



(4.42) 



(4.43) 
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We used Kepler's 3rd lawf to replace R with powers of the orbital frequency O and 
the total mass M — m\ + m 2 . For the purpose of our numerical estimate, we will take 
the members of the binary to have equal masses, so that fi — M/A: 

M 5/3 2/3 

(4.44) 



r 

Finally, we insert numbers corresponding to plausible sources: 

M \ 5 / 3 / 1 hour\ 2 / 3 / 1 kiloparsec\ 



lO" 21 



2M Q J V P J V r 



n 22 f M \ 5/3 / 0.01 sec\ 2/3 /100Megaparsecs\ , 4A . 

= 10 \mM-J {—) { r ) ■ ( " 5) 

The first line corresponds roughly to the mass, distance and orbital period (P = 2ir/tt) 
expected for the many close binary white dwarf systems in our galaxy. Such binaries 
are so common that they are likely to be a confusion limited source of GWs for space- 
based detectors, acting in some cases as an effective source of noise. The second line 
contains typical parameter values for binary neutron stars that are on the verge of 
spiralling together and merging. Such waves are targets for the ground-based detectors 
that have recently begun operations. The tiny magnitude of these waves illustrates 
why detecting GWs is so difficult. 



5. Linearized theory of gravitational waves in a curved background 

At the most fundamental level, GWs can only be defined within the context of an 
approximation in which the wavelength of the waves is much smaller than lengthscales 
characterizing the background spacetime in which the waves propagate. In this section, 
we discuss perturbation theory of curved spacetimes, describe the approximation in 
which GWs can be defined, and derive the effective stress tensor which describes the 
energy content of GWs. The material in this section draws on the treatments given in 
Chapter 35 of Misner, Thorne and Wheeler |J, Sec. 7.5 of Wald |51|. and the review 
articles |3HI32] . 



5.1. Perturbation theory of curved vacuum spacetimes 

Throughout this section we will for simplicity restrict attention to vacuum spacetime 
regions. We consider a one-parameter family of solutions of the vacuum Einstein 
equation, parameterized by e, of the form 

9ab = 9ab + ehab + £ 2 jab + 0(e 3 ) . (5.1) 

Here g^ b is the background metric; it was taken to be the Minkowski metric in Sees. [3 
l4land l2~2"l Here we allow g® b to be any vacuum solution of the Einstein equations. The 
quantity h a b is the linear order metric perturbation, as in the previous sections; j a b is 
a second order metric perturbation which will be used in Sec. 15.31 We can regard e as 
a formal expansion parameter; we set its value to unity at the end of our calculations. 

The derivation of the linearized Einstein equation proceeds as before. Most of the 
formulae for linearized perturbations of Minkowski spacetime continue to apply, with 
r] a b replaced by g^ b , and with partial derivatives d a replaced by covariant derivatives 

fin units with G = 1, and for circular orbits of radius R, R?Q 2 = M. 
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with respect to the background, V B . Some of the formulae acquire extra terms 
involving coupling to the background Riemann tensor. 

Inserting Eq. I|5.1|l into the formula for connection coefficients gives 

n c = \d ad (dc9db + d b9dc - d d g bc ) (5.2) 

= \ (g B ad - eh ad ) (d c g B b + sd c h db + d b g B c 

+ed b h dc - d d g B c - ed d h bc ) + 0{e 2 ) 
= r Ba bc + £dT a bc + 0(e 2 ) . (5.3) 

Here T B bc are the connection coefficients of the background metric g B b , and the first 
order corrections to the connection coefficients are given by 

sr a bc = - \h ad gl^t c + \ g Bad (d c h db + d b h dc - d d h bc ) 

= \g B ad iy B h db + vfh dc - v B h bc ) , (5.4) 

where V B is the covariant derivative operator associated with the background metric. 
Equation (|5.4(l can be derived more easily, at any given point in spacetime, by 
evaluating the expression l|5.2(l in a coordinate system in which the background 
connection coefficients vanish at that point, so that d a = V B . The result (| 5 . 4|) for 
general coordinate systems then follows from general covariance. 

Next, insert the expansion l|5.3(l of the connection coefficients into the formula 

R a bcd — 9 c ^ bd — d d r bc + r" e rj; d — T de ri c (5.5) 

for the Riemann tensor. Evaluating the result in a coordinate system in which 
T B bc = at the point of evaluation gives 

R\ cd = d c T Ba bd ~ d d T Ba bc + e {d c 5T a bd ~ d d STl) + 0(e 2 ) 

= R Ba bcd + eSR a bcd + 0(e 2 ) . (5.6) 

Here R Ba bcd is the Riemann tensor of the background metric, and SR a bcd — d c ST bd — 
d d ST bc is the linear perturbation to the Riemann tensor. It follows from general 
covariance that the expression for SR a bcd in a general coordinate system is 

8R\ cd = V B 5T a bd V B ST a bc . (5.7) 

Using the expression (|5.4I) now gives 

SR a bcd = I h B vfh a d + V B V B ^ fc - V B V Ba h bd 



2 



V?V B ^ C - V B V B c h\ + V B V B a h bc ) . (5.8) 



SR bd = --U B h bd - -V B V B h + W B Wf b h a d) . (5.9) 



Contracting on the indices a and c yields the linearized Ricci tensor 5R b 

where D B = V B V Ba , indices are raised and lowered with the background metric, and 
h = h a a . Reversing the trace to obtain the linearized Einstein tensor 5G bd , and writing 
the result in terms of the trace-reversed metric perturbation 

Kb = h ab -\g B b g Bcd h cd (5.10) 
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yields the linearized vacuum Einstein equation 

= 6G bd = - ±n B h bd + R B dbc h ac 

- \af d ^ B c h ac + V B a h a d + \^ B a h% • (5.H) 

As in Sec. [3] the linearized Einstein equation can be simplified considerably by a 
suitable choice of gauge. Under a gauge transformation parameterized by the vector 
field £ a , the metric transforms as 

h ab -> h' ab = h ab - 2Vf a &) ; (5.12) 
the divergence of the trace-reversed metric perturbation thus transforms as 

V Ba h' ab = V Ba h ab -0 B Z b , (5.13) 
We can enforce in the new gauge the transverse condition 

V B X b -0 (5.14) 

by requiring that satisfy the wave equation B £ b = V Ba h ab . We can further 
specialize the gauge to satisfy h' — 0. Dropping the primes, the metric perturbation 
is thus traceless and transverse: 

V Ba h ab = h = 0. (5.15) 

In this gauge the linearized Einstein equation (|5.11() simplifies to 

= SG bd = -^O B h bd + R B dbc h ac . (5.16) 

(Note, however, that one cannot in this context impose the additional gauge conditions 
ho a = used in the definition of TT gauge for perturbations of flat spacetime.) 

To see that the traceless condition h = can be achieved, note that the trace 
transforms as 

fc-> ti = h- 2V Ba £ Q . (5.17) 

Therefore it is sufficient to find a vector field £ a that satisfies n B £a = and 

V Sa £ a - h/2 = 0. (5.18) 

We can choose initial data for £ on any Cauchy hypersurface for which the quantity 
H5.18fl and also its normal derivative vanish. Since the quantity (|5.18() satisfies the 
homogeneous wave equation by Eqs. (|5.11|) and l|5.14|l . it will vanish everywhere. 

The wave equation l|5.1(j|l differs from its flat spacetime counterpart (|2.1t)|) in 
two respects: First, there is an explicit coupling to the background Riemann tensor; 
and second, there is a coupling to the background curvature through the connection 
coefficients that appear in the covariant wave operator n B . In the limit (discussed 
below) where the wavelength of the waves is much smaller than the lengthscales 
characterizing the background metric, these couplings have the effect of causing 
gradual evolution in the properties of the wave. These gradual changes can be 
described using the formalism of geometric optics, which shows that GWs travel along 
null geodesies with slowly evolving amplitudes and polarizations. See Ref. for 
a detailed description of this formalism. Outside of the geometric optics limit the 
curvature couplings in Eq. (|5.1fci|) can cause the dynamics of the metric perturbation 
to be strongly coupled to the dynamics of the background spacetime. An example of 
such coupling is the parametric amplification of metric perturbations during inflation 
in the early Universe |52j . 
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5. 2. General definition of gravitational waves: The geometric optics regime 

The linear perturbation formalism described in the last section can be applied to 
any perturbations of any vacuum background spacetime. Its starting point is the 
separation of the spacetime metric into a background piece plus a perturbation. 
In most circumstances this separation is merely a mathematical device and can be 
chosen arbitrarily; no unique separation is determined by local physical measurements. 
[Although g^ b and h ab are uniquely determined once one specifies the one parameter 
family of metrics g ab {e), a given physical situation will be described by a single 
metric g a b(£o) for some fixed value of Eq of e, not by the one parameter family of 
metrics.] However, in special circumstances, a unique separation into background 
plus perturbation is determined by local physical measurements, and it is only in this 
context that GWs can be defined. Such circumstances arise when the wavelength A of 
the waves is very much smaller than the characteristic lengthscales £ characterizing 
the background curvature. In this case one can define the background metric and 
perturbation, to linear order, via 

9a b =(g ab ), (5.19) 
eh ab = g ab - g® b . (5.20) 

Here the angular brackets (...) denote an average over lengthscales large compared to 
A but small compared to £; a suitable covariant definition of such averaging has been 
given by Brill and Hartle [121 • A useful analogy to consider is the surface of an orange, 
which contains curvatures on two different lengthscales: An overall, roughly spherical 
background curvature (analogous to the background metric), and a dimpled texture 
on small scales (analogous to the GW). The regime A <C C is called the geometric 
optics regime. 

We will argue below that the short-wavelength perturbation sh ab gives rise to 
an effective stress tensor of order e 2 h 2 /X 2 , where h is a typical size of h ab . This 
effective stress tensor contributes to the curvature of the background metric g^. This 
contribution to the curvature is ^ l/C 2 . It follows that e 2 h 2 /X 2 ^ l/£ 2 , or 

eh<^. (5.21) 

Since we are assuming that A <C £, it follows that the short-wavelength piece eh ab 
of the metric is small compared to the background metric, and so we can use the 
perturbation formalism of Sec. 15. II Consider now the splitting of the Riemann tensor 
into a background piece plus a perturbation given by Eq. I|5.6[) : 

Rabcd = Rabcd + £ $R 

abed + 0(e 2 ). (5.22) 

By the definition H5.19|l of the background metric, it follows that g^ and R^ bcd vary 
only over lengthscales <; £, and therefore it follows that to a good approximation 

(Rabcd) — Rabcd ■ (5.23) 

Hence the perturbation to the Riemann tensor can be obtained via 

zSRabcd = Rabcd ~ (Rabcd) , (5-24) 

the same unique and local procedure as for the metric perturbation H5.20|) . This 
Riemann tensor perturbation is often called the GW Riemann tensor, it is a tensor 
characterizing the GWs that propagates in the background metric g® b . 

The operational meaning of the GW fields eh ab and eSR a bcd follows directly from 
the equivalence principle and from their meaning in the context of flat spacetime (Sec. 
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[21 • Specifically, suppose that V is a point in spacetime, and pick a coordinate system 
in which g\\ = Tj a b and r B ^ c = at V . Then we have 



where x is the distance from V . Therefore, within a spacetime region around V 
in which x <C C, the flat-spacetime perturbation theory and measurement analysis of 
Sec.[3]can be applied. Thus, the gravitational waveforms seen by observers performing 
local experiments will just be given by components of the GW Riemann tensor in the 
observer's local proper reference frames. 

We remark that the splitting of the metric into a background plus a linear 
perturbation can sometimes be uniquely defined even in the regime A ~ C. Some 
examples are when the background spacetime is static (eg perturbations of a static 
star), or homogeneous (eg Friedman- Robertson- Walker cosmological models). In these 
cases the dynamic metric perturbation are not actually GWs, although their evolution 
can be computed using the linearized Einstein equation. For example, consider the 
evolution of a metric perturbation mode which is parametrically amplified during 
inflation in the early Universe. At early times during inflation, the mode's wavelength 
A is smaller than the Hubble scale (£); the mode is said to be "inside the horizon" . Any 
excitation of the mode is locally measurable (although such modes are usually assumed 
to start in their vacuum state). As inflation proceeds, the mode's wavelength redshifts 
and becomes larger due to the the rapid expansion of the Universe, and eventually 
becomes larger than the Hubble scale £; the mode is then "outside the horizon". At 
this point, excitations in the mode are not locally measurable and are thus not GWs. 
Finally, after inflation ends the mode "re-enters the horizon" and excitations of the 
mode are locally measurable. The mode is now a true GW once again. 

Finally, we note that for perturbations of flat spacetime, the definition of GWs 
given here does not always coincide with the definition in terms of the TT component 
of the metric given in Sec. 12.21 However, far from sources of GWs (the regime relevant 
to observations), the two definitions do coincide. This is because the TT piece of the 
metric will vary on scales of a wavelength A which is short compared to the lengthscale 
~ r over which other pieces of the metric vary (except for other dynamic pieces of the 
metric such as the time- varying quadrupolc term in the gauge-invariant field $; those 
pieces vary on short lengthscales but are unimportant since they are smaller than the 
TT piece by a factor ~ A 2 /r 2 or smaller). 

5.3. Effective stress-energy tensor of gravitational waves 

Two major conceptual building blocks are needed for the derivation of the energy and 
momentum carried by GWs The perturbation theory of Sec. 15.11 generalized to 
second order in e, and the separation of lengthscales A <C C discussed in the previous 
subsection. 

We start by discussing the second order perturbation theory. By inserting the 
expansion Ij5.1|l into the vacuum Einstein equation we obtain 




(5.25) 



— G a b 



G ab [gf d ] + eG% [h cd ; gf f ] + e 2 G$ [j cd ; gf f ] + s 2 G% [h cd ; g, 
+ Q(e 3 ) . 



5 



(5.26) 
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Here G ab [gf d ] is the Einstein tensor of the background metric, and G^[...;gf^\ is the 
linear differential operator on metric perturbations giving the linear perturbation to 
the Einstein tensor generated by a metric perturbation. The explicit expression for 
Gab [hcdi gff] is given by Eq. i|5.11|) . The term G^ [h cd ; gf^] is the piece of the Einstein 
tensor that is quadratic in h ab ; it can be computed by extending the computation of 
Sec. 1(5.1(1 to one higher order, and is a sum of terms of the form h ab V^V^h e j and 
( V a^fec)(Vf h ef ) with various index contractions; see Eq. (35.58b) of MTW @|. It's 
worth recalling that j ab is a second order metric perturbation. Wc must take the 
calculation to second order to compute the effective stress energy tensor of the waves 
since an averaging is involved — the first order contribution vanishes by the oscillatory 
nature of the waves. 

Equating to zero the coefficients of the different powers of e we obtain the vacuum 
Einstein equation for the background spacetime 

Gab [fled] = 0, (5.27) 
the linearized Einstein equation 

G a 1 b ) [h cd ;gf f ]=0, (5.28) 
together with the equation for the second-order metric perturbation j ab 

G^ Ha; gf f ] = [h cd ; gf f ] ■ (5.29) 

We now specialize to the geometric optics regime A <C C. We split the second 
order metric perturbation into a piece (j a b) that is slowly varying, and a piece 

Aj ab = jab - (jab) (5.30) 

that is rapidly varying. The full metric can be now be written 

gab = {gab + e 2 (j ab )) + (eh ab + £ 2 Aj ab ) + 0{e 3 ) , (5.31) 

where the first term varies slowly on lengthscales ~ £, and the second term varies 
rapidly on lengthscales ~ A. Consider next the average of the second-order Einstein 
equation ((5.29(1 . Using the fact that the averaging operation (...) commutes with 
derivatives we get 

G^Kjc^gff] = -{G^ih^g*}) . (5.32) 
Subtracting Eq. 1)5.32(1 from Eq. ((5.29(1 gives an equation for Aj a i,: 

G§ [Aj cd ] = -G% [h cd ; g*] + (G$ [h cd ; g* f ]) . (5.33) 
Equation 1(5.32(1 can be rewritten using Eq. 1(5.27(1 asf 

G ab [gf d + s 2 ( Jcd )} = SttT™ + 0(e 3 ) , (5.34) 
where the effective GW stress-energy tensor is 



T ab = --^(G ab J [hcd;g ef }) ■ (5.35) 



In the effective Einstein equation ((5.34(1 , all the quantities vary slowly, on lengthscales 
~ C The left-hand-side is the Einstein tensor of the slowly varying piece of the metric. 
The right hand side is the effective stress energy tensor, obtained by taking an average 

fOur derivation of the effective Einstein equation 15.341 requires the assumption £ 2 (j a b) *K 9^ b , 
since we use second order perturbation theory. However the final result is valid without this 
assumption 1541 : the curvature generated by the GWs can be comparable to the background curvature. 
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of the quadratic piece of the second order Einstein tensor. It follows from Eq. 
that T , ' e is conserved with respect to the metric g^ b + £ 2 (j a b)- I n particular, to 
leading order in e it is conserved with respect to the background metric g^ b . 

The effect of the GWs is thus to give rise to a correction (j ab ) to the background 
metric. This correction is locally of the same order as Aj a f,, the rapidly varying piece 
of the second order metric perturbation. However, any measurements that probe 
only the long-lengthscale structure of the metric (for example measurements of the 
gravitating mass of a radiating source over timescales long compared to A) are sensitive 
only to (jab)- Thus, when one restricts attention to long lengthscales, GWs can thus 
be treated as any other form of matter source in general relativity. Typically (j ab } 
grows secularly with time, while Aj ab does not. 

A fairly simple expression for the effective stress-energy tensor can be obtained 
as follows. Schematically, the effective stress-energy tensor has the form 

TT CS ~ Hh ab V^h ef ) + &((V% c )(V^ e/ )) , (5.36) 

where & means "a sum of terms obtained by taking various contractions of" . In 
this expression gradients scale as 1/A, so VjfV^ ~ 1/A 2 . However the commutator 
of two derivatives scales as the background Riemann tensor, which is of order 1/C 2 . 
Therefore, up to corrections of order A 2 /£ 2 which can be neglected, one can freely 
commute covariant derivatives in the expression 

(jESEJ- Also, the average of any total 
derivative will vanish in the limit A <C C if the averaging lengthscale is taken to be 
V \C. Therefore one can flip derivatives from one factor to another inside the averages 
in Eq. (|5.3(il) . as in integration by parts. Using these manipulations the expression for 
the effective stress-energy tensor simplifies to 0] 02] 



In gauges satisfying the transverse-traceless conditions (|5.15() this reduces to 

T ^ W ' Cff = i ^ h ^ hCd ) ■ ( 5 - 38 ) 
For example, for the plane wave propagating in the z direction in flat spacetime, given 

by 

h X x = — h yy = ho cos(ujt — ojz) , 

h a b = (all other components) , (5.39) 
the energy density and energy flux are given by 

T » = T u = 3f {cos 2 {ujt wz)> = (5 . 40) 

1D7T 32-7T 

If we restore factors of G and c, and insert numbers typical of bursts of waves that we 
hope to detect, we get the energy flux 

T tz = 1.5 mW m" 2 ( " ( ' ( 5 .41) 

V 10- 22 J V lkHz / 

where / = oj/(2it). Note that this is a large energy flux by astronomical standards, 
despite the tiny value of ho', it is comparable to the flux of reflected sunlight from a 
full moon |32| . 
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6. A brief survey of gravitational wave astronomy 

Having now reviewed the basic theory and properties of GWs, we conclude this article 
by very briefly surveying the properties of important potential sources of GWs. Our 
goal is to give some indication of the value that GWs may provide for astronomical 
observations; much of this material is updated from a previous survey article, Ref. 
|29|. We note that since the focus of this article is intended to be the theory GW 
sources (and that this article is significantly longer than was intended or requested), 
we are quite a bit more schematic in our treatment here than we have been in the rest 
of this paper. This final section is intended to be a very brief, somewhat superficial 
survey, rather than a detailed review. 

We begin by contrasting gravitational radiation with electromagnetic radiation, 
which forms the basis for almost all current astronomical observations: 

Electromagnetic waves interact strongly with matter; GWs do not. The weak 
interaction of GWs is both blessing and curse: It means that they propagate from 
emission to Earth-bound observers with essentially zero absorption, making it possible 
to probe astrophysics that is hidden or dark to electromagnetic observations — e.g., 
the coalescence and merger of black holes, the collapse of a stellar core, the dynamics 
of the early Universe. It also means that detecting GWs is very difficult. Also, because 
many of the best sources are hidden or dark, they are very poorly understood today 
— we know very little about what are likely to be some of the most important sources 
of GWs. 

Electromagnetic radiation typically has a wavelength smaller than the size of the 
emitting system, and so can be used to form an image of the source. This is because 
electromagnetic radiation is usually generated by moving charges in the environment 
of some larger source — e.g., an atomic transition in interstellar gas, or emission 
from hot plasma in a stellar environment. By contrast, the wavelength of gravitational 
radiation is typically comparable to or larger than the size of the radiating source. 
GWs are generated by the bulk dynamics of the source itself — e.g., the motion 
of neutron stars in a binary. As a consequence, GWs cannot be used to form an 
image: The radiation simply does not resolve the generating system. Instead, GWs 
are best thought of as analogous to sound — the two polarizations carry a stereophonic 
description of the source's dynamics. 

Gravitons in a gravitational-wave burst are phase coherent; photons in electromagnetic 
signals are usually phase-incoherent. This arises from the fact that each graviton is 
generated from the same bulk motion of matter or of spacetime curvature, while each 
photon is normally generated by different, independent events involving atoms, ions 
or electrons. Thus GWs are similar to laser light. We can take advantage of the 
phase coherence of GWs to enhance their detectability. Matched filtering techniques 
for detecting GW bursts with well-modeled functional form (like those generated by 
coalescing compact binaries) extend the distance to which sources can be seen by a 
factor of roughly the square root of the number of cycles in the waveform, a significant 
gain |35). 

An extremely important consequence of this coherency is that the direct observable 
of gravitational radiation is the strain h, a quantity that falls off with distance as 
1/r. Most electromagnetic observables are some kind of energy flux, and so fall off 
with a 1/r 2 law; measuring coherent GWs is analogous to measuring a coherent, 1/r 
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electromagnetic radiation field. This comparatively slow fall off with radius means 
that relatively small improvements in the sensitivity of GW detectors can have a large 
impact on their science: Doubling the sensitivity of a detector doubles the distance to 
which sources can be detected, increasing the volume of the Universe to which sources 
are measurable by a factor of 8. Every factor of two improvement in the sensitivity 
of a GW observatory should increase the number of observable sources by about an 
order of magnitude. 

In most cases, electromagnetic astronomy is based on deep imaging of small fields of 
view: Observers obtain a large amount of information about sources on a small piece 
of the sky. GW astronomy will be a nearly all-sky affair: GW detectors have nearly 4ir 
steradian sensitivity to events over the sky. A consequence of this is that their ability 
to localize a source on the sky is not good by usual astronomical standards; but, it 
means that any source on the sky will be detectable, not just sources towards which the 
detector is "pointed" . The contrast between the all-sky sensitivity but poor angular 
resolution of GW observatories, and the pointed, high angular resolution of telescopes 
is very similar to the angular resolution contrast of hearing and sight, strengthening 
the useful analogy of GWs with sound. 

From these general considerations, we turn now to specifics. It is useful to 
categorize GW sources (and the methods for detecting their waves) by the frequency 
band in which they radiate. Broadly speaking, we may break the GW spectrum into 
four rather different bands: the ultra low frequency band, 10 _18 Hz <, f <, 10~ 13 Hz; 
the very low frequency band, 10~ 9 Hz <> f <> 10~ 7 Hz; the low frequency band, 
10 -5 Hz <j f <, 1 Hz; and the high frequency band, 1 Hz f <, 10 Hz. 

For compact sources, the GW frequency band is typically related to the source's 
size R and mass M. Here the source size R means the lengthscale over which the 
source's dynamics vary; for example, it could be the actual size of a particular body, 
or the separation of members of a binary. The "natural" GW frequency of such a 
source is few ~ (l/27r)y / GM/i? 3 . Because R ^ 2GM/c 2 (the Schwarzschild radius 
of a mass M), we can estimate an upper bound for the frequency of a compact source: 



This is a rather hard upper limit, since many interesting sources are quite a bit larger 
than 2GM/c 2 , or else evolve through a range of sizes before terminating their emission 
at R ~ 2GM/c 2 . Nonetheless, this frequency gives some sense of the types of compact 
sources that are likely to be important in each band — for example, high frequency 
compact sources are of stellar mass (several solar masses); low frequency compact 
sources are thousands to millions of solar masses, or else contain widely separated 
stellar mass bodies. 

6.1. High frequency 

The high frequency band, 1 Hz <, f <> 10 Hz, is targeted by the new generation 
of ground-based laser interferometric detectors such as LIGO. (It also corresponds 
roughly to the audio band of the human ear: When converted to sound, LIGO 
sources are audible to humans.) The low frequency end of this band is set by the 
fact that it is extremely difficult to prevent mechanical coupling of the detector to 
ground vibrations at low frequencies, and probably impossible to prevent gravitational 




(6.1) 
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coupling to ground vibrations, human activity, and atmospheric motions |46l 1471 l4"5| . 
The high end of the band is set by the fact that it is unlikely any interesting GW 
source radiates at frequencies higher than a few kilohertz. Such a source would have 
to be relatively low mass ( ^ lM©) but extremely compact [cf. Eq. Q6.ip ]. There 
are no known theoretical or observational indications that gravitationally collapsed 
objects in this mass range exist. 

The article by Danzmann in this volume discusses the detectors relevant to 
this frequency band in some detail; our discussion here is limited to a brief survey 
of these instruments. Several interferometric GW observatories are either operating 
or being completed in the United States, Europe, Japan, and Australia. Having 
multiple observatories widely scattered over the globe is extremely important: The 
multiplicity gives rise to cross-checks that increase detection confidence and also aids 
in the interpretation of measurements. For example, sky location determination and 
concomitant measurement of the distance to a source follows from triangulation of 
time-of-flight differences between separated detectors. 

The major interferometer projects are: 

• LIGO. The Laser Interferometer Gravitational-wave ObservatoryjHE] consists 
of three operating interferometers: A single four kilometer interferometer in 
Livingston, Louisiana, as well as a pair of interferometers (four kilometers and two 
kilometers) in the LIGO facility at Hanford, Washington. The sites are separated 
by roughly 3000 kilometers, and are situated to support coincidence analysis of 
events. 

• Virgo. Virgo is a three kilometer French-Italian detector under construction 
near Pisa, Italy In most respects, Virgo is quite similar to LIGO. A major 
difference is that Virgo employs a very sophisticated seismic isolation system that 
promises extremely good low frequency sensitivity. 

• GEO600. GEO600 is a six hundred meter interferometer constructed by a 
German-English collaboration near Hanover, Germany Despite its shorter 
arms, GEO600 achieves sensitivity comparable to the multi-kilometer instruments 
using advanced interferometry techniques. This makes it an invaluable testbed 
for technology to be used in later generations of the larger instruments, as well 
as enabling it to make astrophysically interesting measurements. 

• TAMA300. TAMA300 is a three hundred meter interferometer operating near 
Tokyo. It has been in operation for several years now [SHj. The TAMA team 
is currently designing a three kilometer interferometer |59| . building on their 
experiences with the three hundred meter instrument. 

• ACIGA. The Australian Consortium for Interferometric Gravitational- Wave 
Astronomy is currently constructing an eighty meter research interferometer near 
Perth, Australia (HOj, hoping that it will be possible to extend it to multi-kilometer 
scale in the future. Such a detector would likely be a particularly valuable addition 
to the worldwide stable of detectors, since all the Northern Hemisphere detectors 
lie very nearly on a common plane. An Australian detector would be far outside 
this plane, allowing it to play an important role in determining the location of 
sources on the sky. 

The LIGO, GEO, and TAMA instruments have now been operating for several 
years; see Refs. ^| El El 120] f° r the results and upper limits from the first set of 
observations. All of these detectors have or will have sensitivities similar to that 
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Figure 2. Sensitivity goals of the initial LIGO interferometers, and facility limits 
on the LIGO sensitivity (taken from Ref. I25I V 



illustrated in Fig. [5] (which shows, in particular, the sensitivity goal of the first 
generation of LIGO interferometers). This figure also shows the "facility limits" - 
the lowest noise levels that can be achieved even in principle within an interferometer 
facility. The low level facility limits come from gravity- gradient noise: noise arising 
from gravitational coupling to fluctuations in the local mass distribution (such as 
from seismic motions in the earth near the test masses |46| . human activity near the 
detector |47|. and density fluctuations in the atmosphere |48|1. At higher frequencies, 
the facility limit arises from residual gas (mostly hydrogen) in the interferometer 
vacuum system — stray molecules of gas effectively cause stochastic fluctuations in 
the index of refraction. 

We now survey the more well-understood possible sources of measurable GWs in 
the high-frequency band. We emphasize at this point that such a listing of sources 
can in no way be considered comprehensive: We are hopeful that some GW sources 
may surprise us, as has been the case whenever we have studied the Universe with a 
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new type of radiation. 



6.1.1. Coalescing compact binaries Compact binaries — binary star systems in which 
each member is a neutron star or black hole — are currently the best understood 
sources of GWs. Double neutron stars have been studied observationally since the 
mid 1970s; five such systems ^1 El [HI El E! tight enough to merge within a few 
10 8 or 10 9 years have been identified in our Galaxy. Extrapolation from these observed 
binaries in the Milky Way to the Universe at large |^ 1621 1631 R^f] indicates that GW 
detectors should measure at least several and at most several hundred binary neutron 
star mergers each year (following detector upgrades; the expected rate for initial 
detectors is of order one event per several years, so that measurement of an event 
is plausible but of fairly low probability) . Population synthesis (modeling evolution of 
stellar populations) indicates that the measured rate of binaries containing black holes 
should likewise be interestingly large (perhaps even for initial detectors) 65 66, 67JEHI • 
The uncertainties of population synthesis calculations are rather large, however, due 
to poorly understood aspects of stellar evolution and compact binary formation; data 
from GW detectors is likely to have a large impact on this field. 



6.1.2. Stellar core collapse Core collapse in massive stars (the engine of Type II 
supernova explosions) has long been regarded as likely to be an important source 
of GWs; see, for example, Ref. jEH for an early review. Stellar collapse certainly 
exhibits all of the necessary conditions for strong GW generation — large amounts of 
mass (1 — 100 M Q ) flow in a compact region (hundreds to thousands of kilometers) 
at relativistic speeds (v/c <; 1/5). However, these conditions are not sufficient to 
guarantee strong emission. In particular, the degree of asymmetry in collapse is not 
particularly well understood. 

If the core of a star is very rapidly rotating during collapse, then instabilities 
may develop which lead to strong GW emission [70]. If such instabilities develop, core 
collapse GWs could be detected from events as far away as 10 Megaparsecs [TT], a 
distance encompassing enough galaxies that several events per year would be likely. 
Most models of massive stars, however, indicate that such rapid rotation is not likely 
(e.g., |7_3). Even without the growth of instabilities, the asymmetric dynamics of core 
collapse is likely to lead to wave emission which would be detectable within the Local 
Group of galaxies, with perhaps an event every few years detectable by advanced 
interferometers [73j . The wave strength is likely to correlate strongly with the degree 
of asymmetry in the supernova. If the GW event has an electromagnetic or neutrino 
counterpart we may gain a wealth of knowledge regarding the state of the precollapse 
core |71) . 



6.1.3. Periodic emitters Periodic sources of GWs radiate at constant or nearly 
constant frequency, like radio pulsars. In fact, the prototypical source of continuous 
GW is a rotating neutron star, or GW pulsar. A non-axisymmetry in a neutron star 
crust (caused, for example, by an oblateness that is misaligned with the star's spin 
axis) will radiate GWs with characteristic amplitude 

h --, J r- (6-2) 

Here I is the star's moment of inertia, / is the wave frequency, r is the distance to the 
source, and e is the dimensionless fractional distortion e = (I xx — I yy )/I, where Jy 
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is the moment of inertia tensor. The crucial parameter e characterizes the degree to 
which the star is distorted; it is rather poorly understood. Upgraded interferometers in 
LIGO could set an upper limit on e of order 10 -6 for sources at ~ 10 kpc fZ7\. Various 
mechanisms have been proposed to explain how a neutron star can be distorted to 
give a value of e that is interesting as a GW source; see [23 ESI for further discussion. 
Examples of some interesting mechanisms include misalignment of a star's internal 
magnetic field with the rotation axis [77) and distortion by accreting material from a 
companion star [JHUZHj (discussed in more detail below). 

Whatever the mechanism generating the distortion, it is clear that e will be small, 
so that h ~ 10 -24 or smaller — quite weak. Measuring these waves will require 
coherently tracking their signal for a large number of wave cycles. Coherently tracking 
N cycles boosts the signal to noise ratio by a factor ~ y/~N. This is actually fairly 
difficult, since the signal is strongly modulated by the Earth's rotation and orbital 
motion, "smearing" the waves' power across multiple frequency bands. Searching for 
periodic GWs means demodulating the motion of the detector, a computationally 
intensive problem since the modulation is different for every sky position. Unless 
one knows in advance the position of the source, one needs to search over a huge 
number of sky position "error boxes", perhaps as many as 10 13 . One rapidly becomes 
computationally limitedf. (Radio pulsar searches face this same problem, with the 
additional complication that radio pulses are dispersed by the interstellar medium. 
However, radio observers usually use shorter integration times, and often target their 
searches on small regions of the sky, so their computational cost is usually not as 
great.) For further discussion, see [HO]; f° r ideas about doing hierarchical searches 
that require less computer power, see [81) . 

As mentioned above, one particularly interesting mechanism for distorting a 
neutron star comes from accretion of material from a companion star. Accretion 
provides a spin-up torque to a neutron star, 

(dJ/dt) spin _ ap ~ R 2 n*M (6.3) 

(where J is the spin angular momentum, fi* is the orbital frequency of the accreting 
matter as it plunges onto the star, R is the star's radius, and M is the mass 
accretion rate). Without any kind of braking mechanism, the neutron star would 
presumably spin-up until it reaches the "breakup limit", i.e., the spin frequency at 
which centrifugal forces would begin to break it apart; the breakup frequency is 
typically around 2000 - 3000 Hz. 

Observations have shown |82) that accreting neutron stars do, in fact, appear to 
have a "speed limit" — no accreting neutron star has been observed to spin faster 
than 619 Hz [53]. This is consistent with the fact that the fastest radio pulsar! has 
a spin period of 641 Hz [84|. This suggest that some mechanism is removing angular 
momentum from the neutron star. A plausible and very attractive possibility for how 

fThis rather large number of patches on the sky is driven by the possible need to search for 
high frequency pulsars over several months of observation. The difference A/ between the Doppler 
frequency shifts for two adjacant sky patches separated by an angle 59 is of order A/ ~ v^fSB/c, 
where ~ 3 X 10 4 ms~ 1 is the Earth's orbital velocity and / is the gravitational wave frequency. 
The phase error over an observation time T t, s is of order A/T t, s . Demanding that this be less 
than unity yields 69 < c/(v^fT). The number of independent sky patches is then N p ~ 4,nS9~ 2 ~ 
4tto 2 / 2 T o 2 bs /c 2 ~ 10 13 for / = 1000 Hz and T obs = 1/3 year. Fewer positions would be needed 
if either the maximum frequency or the integration time is reduced; the figures given here set the 
maximum values that are plausible. See Ref. 1801 for more details. 

JThe so-called "recycled" radio pulsars spin at frequencies ~ several hundred Hertz; they are 
believed to be the fossils of accreting neutron stars. 
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this angular momentum is removed is via GW emission. Because the spin-down torque 
due to GW emission grows sharply with spin frequency, 

(rf^M) spin-down K fi5 (Quadrupole emission) , (6.4) 

the limiting spin obtained by balancing the torques 1)6. 3JI and l|6. 4JI is relatively 
insensitive to the mass accretion rate M. Such a mechanism was originally suggested 
by Wagoner [53], and was revived by Bildsten [7SJ to explain the narrow clustering 
in spin frequency of accreting low-mass x-ray binaries (LMXBs) . Various mechanisms 
could provide the spin-down torque — Bildsten originally suggested a quadrupole 
moment in the spinning star could be induced by a thermally varying electron capture 
mechanism, but also noted that the r-mode instability (see, e.g., [5S] for a review) could 
be excited, leading to a similar spin-down law. Whatever the mechanism, accreting 
neutron stars are obvious and very attractive targets for observing campaigns with 
GW detectors, particularly given that their sky positions are well known. 



6.1. 4-. Stochastic backgrounds Stochastic backgrounds are "random" GWs, arising 
from a large number of independent, uncorrelated sources that are not individually 
resolvable. A particularly interesting source of stochastic waves is the dynamics of the 
early Universe, which could produce an all-sky GW background, similar to the cosmic 
microwave background; see Refs. [S7J [8J3 EH f° r detailed reviews. Stochastic waves 
can be generated in the early Universe via a variety of mechanisms: amplification of 
primordial fluctuations in the Universe's geometry via inflation, phase transitions as 
previously unified interactions separate, a network of vibrating cosmic strings, or the 
condensation of a brane from a higher dimensional space. These waves can actually 
extend over a wide range of frequency bands; waves from inflation in particular span 
all bands, from ultra low frequency to high frequency. 

Stochastic backgrounds are usually idealized as being stationary, isotropic and 
homogeneous. They are thus characterized by their energy density per unit frequency, 
dp gw /df. This is often parameterized in terms of the energy density per unit 
logarithmic frequency divided by the critical energy density to close the Universe 

1 

p cr it din/ 

where p cr it = 3Hq/8ttG is the critical density and Hq is the value of the Hubble 
constant today. Different cosmological sources produce different levels of f2 gw (/), 
centered in different bands. In the high frequency band, waves produced by inflation 
are likely to be rather weak: estimates suggest that the spectrum will be flat across 
LIGO's band, with magnitude Sl gw ~ 10~ 15 at best [HO]- Waves from phase transitions 
can be significantly stronger, but are typically peaked around a frequency that depends 
on the temperature T of the phase transition |871 IOT] : 

^ k ~ 100Hz (To5w) ' (6 - 6) 

Because of their random nature, stochastic GWs look just like noise. Ground- 
based detectors will measure stochastic backgrounds by comparing data at multiple 
sites and looking for "noise" that is correlated 88, 92 . For comparing to a detector's 
noise, one should construct the characteristic stochastic wave strain, 



M/) = -^> (6-5) 



ft(xr 3 /yfi gw (/)A/, (6.7) 
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where A/ is the frequency band across which the measurement is made. For further 
discussion and the proportionality constants, see Note that if O gw (/) is constant, 
this strain level grows sharply with decreasing frequency — the most interesting limits 
are likely to be set by measurements at low frequencies. 

Early detectors will have fairly poor sensitivity to the background, constraining 
it to a level f2 gw ~ 5 x 1CP 6 in a band from about 100 Hz to 1000 Hz. This is barely 
more sensitive than known limits from cosmic nucleosynthesis |87j . Later upgrades 
will be significantly more sensitive, able to detect waves with fl gw ~ 10~ 10 , which is 
good enough to place interesting limits on backgrounds from some phase transitions. 

6. 2. Low frequency 

There is no hope of measuring GWs in the low frequency band, 10 -5 Hz < / <, 1 Hz, 
using a ground-based instrument: Even if it were possible to completely isolate one's 
instrument from local ground motions, gravitational coupling to fluctuations in the 
local mass distribution ultimately limits the sensitivity to frequencies / <; 1 Hz. 
Nonetheless, many extremely interesting sources radiate in this band. The only way 
to measure these waves is to build a detector in the quiet environment of space, far 
removed from low- frequency noise sources. 

Such an instrument is currently being designed jointly by NASA in the United 
States and ESA, the European Space Agency: LISA, the Laser Interferometer Space 
Antenna. If all goes well, LISA will be launched into orbit in 2013 or so. LISA will 
be a laser interferometer, similar in concept to the ground-based detectors: Changes 
in the distance between widely separated test masses will be monitored for GWs. 
However, LISA's scale is vastly larger than that of its ground-based cousins, and so 
details of its operations are quite different. In particular, LISA will have armlengths 
L ~ 5 x 10 6 km. The three spacecraft which delineate the ends of LISA's arms are 
placed into orbits such that LISA forms a triangular constellation orbiting the sun, 
inclined 60° with respect to the plane of the ecliptic and following the Earth with a 
20° lag. This configuration is sketched in Fig. [3] Since it essentially shares Earth's 
orbit, the constellation orbits the sun once per year, "rolling" as it does so. This 
motion plays an important role in pinpointing the position of sources by modulating 
the measured waveform — the modulation encodes source location and makes position 
determination possible. 

The 3 spacecraft each contain two optical assemblies, each of which houses a 
1 Watt laser and a 30 centimeter telescope. Because of the extreme lengths of the 
interferometer's arms, Fabry-Perot interferometry as in the ground-based detectors is 
not possible: Diffraction spreads the laser beam over a diameter of about 20 km as 
it propagates from one spacecraft to the other. A portion of that 20 km wavefront is 
sampled with the telescope. That light is then interfered with a sample of light from 
the on-board laser. Each spacecraft thus generates two interference data streams; 
six signals are generated by the full LISA constellation. From these six signals, we 
can construct the time variations of LISA's armlengths and then build both GW 
polarizations. More information and details can be found in 93 . 

Note that the LISA armlengths are not constant — as the constellation orbits, 
the distances between the various spacecraft vary by about 1% (including effects such 
as planetary perturbations). These variations are far larger than the displacements 
produced by GWs, which are of order picometers. However, these variations occur 
over timescales of order months, and are extremely smooth and well modeled. It will 
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Figure 3. Orbital configuration of the LISA antenna. 



not be difficult to remove them from the data leaving clean data in the interesting 
frequency band. Picometer scale variations are not too difficult to measure in this 
band by gathering photons for a time 10 sec <, T 1 day. Even though the bulk of 
the laser's emitted power is lost by diffraction, enough photons are gathered on this 
timescale that the phase shift due to the GW can be determined. 

The GW signals are actually read out by monitoring the position of the so- 
called "gravitational sensor" on each optical assembly; in particular, the position of a 
"proof mass" which floats freely and constitutes the test mass for the LISA antenna 
is monitored. Because it is freely floating, the proof mass follows a geodesic of the 
spacetime. Micronewton thrusters keep the bulk spacecraft centered on these proof 
masses, forcing the craft to follow their average trajectory. In this way, the spacecraft 
are isolated from low frequency forces that could impact the ability to measure GWs 
(e.g., variations in solar radiation pressure). 

We now take a quick tour through some interesting LISA sources: 

6.2.1. Periodic emitters In the high-frequency band, the source of most periodic 
GWs is expected to be isolated neutron stars. LISA's periodic GWs will come 
primarily from binary star systems in the Milky Way, primarily close white dwarf 
binaries. Most of these systems do not generate waves strong enough to backreact 
significantly, so that their frequencies do not change measurably over the course 
of LISA observations. Certain systems are well-known in advance to be sources of 
periodic waves for the LISA band. These sources are understood well enough from 
optical observations that they may be regarded as "calibrators" — LISA had better 
detect them, or else something is wrong! 

Aside from these sources that are known in advance, it is expected that LISA will 
discover a good number of binary systems that are too faint to detect with telescopes. 
Joint observations by LISA and other astronomical instruments are likely to be more 
fruitful than observations with a single instrument alone. For example, it is typically 
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difficult for telescopes to determine the inclination of a binary to the line of sight 
(a quantity needed to help pin down the masses of the binary's members). GWs 
measure the inclination angle almost automatically, since this angle determines the 
relative magnitude of the polarizations h+ and h x - 

The total number of periodic binaries radiating in LISA's band is expected 
to be so large that they will constitute a confused, stochastic background at low 
frequencies — there are likely to be several thousand galactic binaries radiating in 
each resolvable frequency bin. This background will constitute a source of "noise" 
(from the standpoint of measuring other astrophysical sources) that is largely than 
that intrinsic to the instrument noise at / < 10~ 3 Hz. 

6.2.2. Coalescing binary systems containing black holes Coalescing binary black hole 
systems will be measurable by LISA to extremely large distances — essentially to the 
edge of the observable Universe. Even if such events are very rare, the observed 
volume is enormous so an interesting event detection rate is very likely. One class 
of such binaries consists of systems in which the member holes are of roughly equal 
mass (~ 10 5 — 10 8 M Q ). These binaries can form following the merger of galaxies 
(or pregalactic structures) containing a black holes in their cores. Depending on the 
mass of the binary, the waves from these coalescences will be detectable to fairly 
large redshifts (z ~ 5 — 10), possibly probing an early epoch in the formation of the 
Universe's structure 

The other major class of binary systems consists of relatively small bodies (black 
holes with mass ~ 10 - 100 Mq, neutron stars, or white dwarfs) that are captured 
by larger black holes (M ~ 10 5 - 10 7 Af Q ). These "extreme mass ratio" binaries 
are created when the smaller body is captured onto an extremely strong field, highly 
relativistic orbit, generating strong GWs. Such systems are measurable to a distance 
of a few Gigaparsecs if the inspiralling body is a 10 Mq black hole, and to a distance 
of a few hundred Megaparsecs if the body is a neutron star or white dwarf. LISA will 
measure the waves that come from the last year or so of the smaller body's inspiral, and 
thence probe the nature of the larger black hole's gravitational field deep within the 
hole's potential. The rates for such events are not too well understood and depend 
on the details of stellar dynamics in the cores of galaxies. Extremely conservative 
estimates typically find that the rate of measurable events for LISA should be at 
least several per year |95l I96| . Recent thinking suggests that these rates are likely 
to be rather underestimated — black holes (which are measurable to much greater 
distances) are likely to dominate the measured rate, perhaps increasing the rate to 
several dozen or several hundred per year |97) . 

Finally, it is worth noting that many events involving intermediate mass black 
holes — those with masses in the band running from a few 10 2 to a few 10 5 Mq - 
would generate GWs in LISA's sensitive band. There is a large body of tentative 
evidence for the existence of black holes in this mass band (see, e.g., |HEJ for a review), 
though as yet we have no "smoking gun" unambiguous signature for such a hole. If 
such black holes do exist and undergo mergers in sufficient numbers, measurement of 
their waves will make possible a wealth of interesting tests of relativity and could 
untangle some of the mysteries surrounding supermassive black hole formation and 
growth. 
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6.2.3. Stochastic backgrounds As discussed in Sec. 16. l.H ground-based detectors can 
measure a stochastic background by correlating the data streams of widely separated 
detectors. LISA will use a slightly different technique: by combining its six data 
streams in an appropriate way, one can construct an observable that is completely 
insensitive to GWs, measuring noise only |100| . This makes it possible to distinguish 
between a noise-like stochastic background and true instrumental noise, and thereby 
to learn about the characteristics of the background . 

The sensitivity of LISA will not be good enough to set interesting limits on an 
inflationary GW background: LISA will only reach f2 gw ~ 10~ n , about four orders 
of magnitude too large to begin to say something about inflation jSU]. However, 
LISA's band is well placed for other possible sources of cosmological backgrounds. In 
particular, an electroweak phase transition at temperature T ~ 100— 1000 GeV would 
generate waves in LISA's band [cf. Eq. (|6.6|) ]. These waves are likely to be detectable 
if the phase transition is strongly hrst order (a scenario that does not occur in the 
standard model, but is conceivable in extensions to the standard model |91p. 

6. 3. Very low frequency 

The very low frequency band, 10 -9 Hz <J / <J 10 -7 Hz, corresponds to waves with 
periods ranging from a few months to a few decades. Our best limits on waves in 
this band come from observations of millisecond pulsars. First suggested by Sazhin 
|1()2) and then carefully analyzed and formulated by Detweiler |l(K-t| . GWs can drive 
oscillations in the arrival times of pulses from a distant pulsar. Millisecond pulsars 
are very good "detectors" for measurements in this band because they are exquisitely 
precise clocks. The range of frequencies encompassed by the very low frequency band 
is set by the properties of these radio pulsar measurements: the high end of the 
frequency band comes from the need to integrate the radio pulsar data for at least 
several months; the low end comes from the fact that we have only been observing 
millisecond pulsars for a few decades. (One cannot observe a periodicity shorter than 
the span of one's dataset.) A recent upper limit derived from this technique is j!04| 

Ogw/iioo < 9-3 x 10" 8 , 4 x 10~ 9 Hz < / < 4 x 10~ 8 Hz (6.8) 

(where the limit is a 95% confidence limit and h W Q is the Hubble constant in units of 
100 km sec" ^pc -1 ). 

The upper limit (|6.8|l already places constraints on some cosmological models (in 
particular those involving cosmic strings) . With further observations and the inclusion 
of additional pulsars in the datasets, it is likely to improve quite soon. It is possible 
that the background in this band will be dominated by many unresolved coalescing 
massive binary black holes |105j — binaries that are either too massive to radiate in 
the LISA band, or else are inspiralling towards the LISA band en route to a final 
merger several centuries or millenia hence. Constraints from pulsar observations in 
this band will remain an extremely important source of data on stochastic waves in 
the future — the limits they can set on fi gw are likely to be better than can be set by 
any of the laser interferometric detectors. 

6.4- Ultra low frequency 

The ultra low frequency band, 10~ 18 Hz <, f <, 10~ 13 Hz, is better described by 
converting from frequency to wavelength: For these waves, lO -5 ^^ 1 < A < -ff^ 1 , 
where Hq 1 ~ 10 10 light years is the Hubble length today. Waves in this band oscillate 
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on scales comparable to the size of the Universe. They are most likely to be generated 
during inflation: Quantum fluctuations in the spacetime metric are parametrically 
amplified during inflation to relatively high amplitude. The rms amplitude to which 
the waves are amplified depends upon the energy scale £"i n fl of inflation: 



where mp is the Planck mass. Measuring these GWs would be a direct probe of 
inflationary physics, and would determine the inflation energy scale, which is currently 
unknown to within many orders of magnitude. These waves have been described as 
the "smoking gun" signature of inflation |106| . 

During inflation, quantum fluctuations impact both the scalar field which drives 
inflation (the inflaton </>) and the metric of spacetime. There exist independent scalar 
fluctuations (coupled fluctuations in the inflaton and scalar-type fluctuations in the 
metric) and tensor fluctuations (tensor- type fluctuations in the metric). The Fourier 
modes of these scalar and tensor perturbations are describable as harmonic oscillators 
in the expanding Universe |107| . Each mode undergoes zero-point oscillations in the 
harmonic potential. However, the potential itself is evolving due to the expansion 
of the Universe. The evolution of this potential parametrically amplifies these zero- 
point oscillations, creating quanta of the field [57] • During inflation, the Universe's 
scale factor a(t) grows faster than the Hubble length H" 1 , and so each mode's 
wavelength likewise grows faster than the Hubble length. The mode's wavelength 
eventually becomes larger than the Hubble length, or the mode "leaves the horizon" . 
After inflation ends the mode subsequently renters the horizon. For gravitational 
perturbations the number of quanta generated in the mode is proportional to the 
factor by which the Universe expands between the two different horizon crossings. 
Fluctuations in the inflaton seed density fluctuations, Sp(r) = 5<fi(r)(dV/d(f>) [where 
V(<p) is the potential that drives the inflaton field]. The tensor-type fluctuations in 
the spacetime metric are GWs. 

Density fluctuations and GWs both leave an imprint upon the cosmic microwave 
background (CMB). First, each contributes to the CMB temperature anisotropy. 
However, even a perfectly measured map of temperature anisotropy cannot really 
determine the contribution of GWs very well because of cosmic variance: Since we 
only have one Universe to observe, we are sharply limited in the number of statistically 
independent influences upon the CMB that we can measure. Large angular scales are 
obviously most strongly affected by this variance, and these scales are the ones on 
which GW most importantly impact the CMB 108 . 

Fortunately, the scalar and tensor contributions also impact the polarization of 
the CMB. These two contributions can be detangled from one another in a model- 
independent fashion. This detangling uses the fact that the polarization tensor 
P a b(n) on the celestial sphere can be decomposed into tensor harmonics. These 
harmonics come in two flavors, distinguished by their parity properties: The "E- 
modes" or "gradient-type" harmonics Yy m ^ ab (ii) [which pick up a factor (—1)' under 
n — * — n], and the "B-modes" or "curl-type" harmonics Y(j m ) ab {n) [which pick up a 

factor (— 1)' +1 under n — > — n]. These harmonics are constructed by taking covariant 
derivatives on the sphere of the "ordinary" spherical harmonics Yj m (n); see |1()9| for 
details. Because scalar perturbations have no handedness, they only induce gradient- 
type polarization. GWs induce both gradient- and curl-type polarization. Thus, an 
unambiguous detection of the curl-type polarization would confirm production of GWs 




(6.9) 
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by inflation. (A caveat is that is that gravitational lensing can convert E-modes to 
B-modes; this so-called "cosmic shear" ultimately limits the sensitivity to GWs of 
CMB polarization studies |11U| -1 

7. Conclusion 

This article has summarized many of the most important topics in the theory of 
GWs. Due to space and time limitations, we sadly were not able to cover all topics 
with which students of this held should be familiar. In particular, we had hoped to 
include a discussion of strong field relativity and GW emission. We confine ourselves, 
in this conclusion, to a (very) brief discussion of important aspects of this subject for 
GW science, as well as pointers to the relevant literature. 

Linearized theory as described in Secs.|21and|Slis entirely adequate to describe the 
propagation of GWs through our Universe and to model the interaction of GWs with 
our detectors. In some cases, it is even adequate to describe the emission of waves 
from a source, as described in Sec. Q] (although for sources with non-negligible self 
gravity such as binary star systems one has to augment linearized theory as described 
in Sec. l4.2[l . However, many sources have very strong self gravity where the linearized 
treatment is completely inadequate. A variety of formalisms have been developed to 
handle these cases. 

• Post-Newtonian (PN) theory. PN theory is one of the most important of 
these formalisms, particularly for modeling binary systems. Roughly speaking, 
PN theory analyzes sources using an iterated expansion in two variables: The 
"gravitational potential", ~ M/r, where M is a mass scale and r characterizes 
the distance from the source; and velocities of internal motion, v. (In linearized 
theory, we assume <p is small but place no constraints on v.) Newtonian gravity 
emerges as the first term in the expansion, and higher order corrections are found 
as the expansion is iterated to ever higher order. Our derivation of the quadrupole 
formula in Sec. l4.2l givcs the leading order term in the PN expansion of the emitted 
radiation. See Luc Blanchet's recent review and references therein for a 
comprehensive introduction to and explication of this subject. 

• Numerical relativity. Numerical relativity seeks to directly integrate Einstein's 
equations on a computer. Ideally, we would like to to use a well-understood 
model of a GW source (e.g., a binary system in which the field strengths are 
small enough that it is well described by post-Newtonian theory) as "initial data" , 
and then numerically evolve the Einstein equations from that point to some final 
equilibrium configuration. The form in which we normally encounter Einstein's 
equation in textbooks is not well suited to this task — the coordinate freedom of 
general relativity means that there is no notion of "time" built into the equation 
G a b = S>nT ao . One must introduce some notion of time for the concept of "initial 
data" to have any meaning. The 4 dimensions of spacetime are then split into 
3+1 dimensions of space and time. Having made this choice, Einstein's equations 
take on a particular form which is amenable to numerical computation. 

A detailed discussion of numerical relativity is given in the contribution by 
Choptuik to this volume |112j : we also recommend the reviews by Lehner |113| 
and by Baumgarte and Shapiro |37|. For the purpose of our present discussion, 
it suffices to remark that it has proven to be extremely difficult to model some of 
the most interesting and important GW sources. In particular, the final stage of 
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binary black hole mergers — regarded by many as the "Holy Grail" of numerical 
relativity — has proven to be quite a challenge. 

• Perturbation theory. In some cases, GW sources can be modeled as nearly, but 
not quite, identical to some exact solution of the Einstein field equations. For 
example, the end state of binary black hole coalescence must be a single black 
hole. As we approach this final state, the system will be well-modeled as the Kerr 
black hole solution, plus some distortion that radiates away. Another example 
is a binary consisting of a stellar mass compact body orbiting a massive black 
hole. The binary's spacetime will be well-described as a single black hole plus a 
perturbation due to the captured body. These cases can be nicely described using 
perturbation theory: We treat the spacetime as some exact background, g^ b , plus 
a perturbation h ab : 

9ab=gab + h ab- (7.1) 

We are in the perturbative regime if ||/ia&||/||<7abll ^ This system can then be 
analyzed by expanding the Einstein equations for this metric and keeping terms 
to first order in h a b (see Sec. I5.1l for details but without the matter source terms 
included) . 

This approach has proven to be particularly fruitful when the background 
spacetime is that of a black hole. For the case of a Schwarzschild background, the 
derivation of the full perturbation equations is rather straightforward; Rezzolla 
gives a particularly compact and readable summary |114j . Perturbations of Kerr 
black holes are not nearly so simple to describe, largely due to the lack of spherical 
symmetry — expanding the metric as in Eq. I|7.1|l does not prove to be so fruitful 
as it is in the Schwarzschild case. Somewhat miraculously, it turns out that 
progress can be made by expanding the curvature tensor: By expanding the 
Riemann tensor as R a bcd — Rabcd + ^Rabcd and taking an additional derivative of 
the Bianchi identity, 

^eRabcd + V ' dRabec + V c R a bde = , (7.2) 

one can derive a wave-like equation for the perturbation 5R a bcd- This analysis 
was originally performed by Teukolsky; see his original analysis |115| for details. 
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Appendix A. Existence of TT gauge in local vacuum regions in linearized 
gravity 

In this appendix we show that one can always find TT gauges in local vacuum regions 
in linearized gravity. More precisely, suppose that V is a connected open spatial region, 
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and (to,tx) is an open interval of time. Then one can find a gauge on the product 
1Z = (t ,ti) x V that satisfies h tt — h ti — S^hij — djhij — 0, as long as T a0 = 
throughout 1Z. 

The proof involves a generalization of the gauge-invariant formalism of Sec. 12.21 
to finite spacetime regions. We define a decomposition of the metric perturbation h a b 
in terms of quantities 0, 7, hj^ , H, £4 and A using the same equations (|2.29|l - 
l|2.35[) as before. However we replace the boundary conditions <|2.3GI) with 



l\dv = / dt<t>\dv, (A- 1 ) 

Jto 

A|av = 0, (A.2) 

(V 2 A)| 9V =H {dv , (A.3) 

(n x e)|av = 2 / eft (n x /% v , (A.4) 

Jto 



where n is the unit outward-pointing unit normal to dV. The reason for this particular 
choice of boundary conditions will be explained below. These boundary conditions 
define a unique decomposition of the metric within 1Z. 

Next, we compute how the variables </>, /3j, 7, H, Si and A transform under 
general gauge transformations. We use the same parameterization lj2.37J) of the gauge 
transformation as before, except we impose now the boundary condition C^y = 0. 
We find that the transformation laws l|2.38|l - 12.44f> are replaced by the following 



equations which contain some extra terms: 

<p -> 4>-A, (A.5) 

A -> fli-Bi- d l ip , (A.6) 

7 -» -y-A-C + il>, (A.7) 

H -> H-2V 2 C, (A.8) 

A -> A-2C*, (A.9) 

-» £< - 25* + 2rji - 2(t - t )^ , (A.10) 

hf? -» ft£ T - 2%^-) +2(t-t )d i d j il>- (ATI) 

Here ^ is the time-independent, harmonic function defined by V 2 ?/; = and V'lav = 



-^|9V,t=t - Similarly 77, is the time- independent, harmonic transverse vector defined by 
V 2 ?7i = and (n x J7)| 9V = n x B| SV ,t=t - 

We define the variables and Sj by the same equations (|2.45(l - (|2.47|l as 
before. From the transformation laws (|A.5|I - (|A.llJl these variables are still gauge 
invariant, while hJ 3 T is no longer gauge invariant in the present context. Next, 
imposing the linearized vacuum Einstein equations using the expressions (|2.59|l - (|2.6ff) 
yields 

v 2 e = 0, v 2 s, = -2(3*6, v 2 $ = (A.12) 

in V. The boundary conditions (|A.lJl - l|A.4|) together with the definitions l|2.45Jl - 
H2.47|l imply that the boundary conditions on the gauge invariant variables are 

$|9V = ®\ov = S|av - 0. (A.13) 

(This is why we choose those particular boundary conditions.) Therefore all the gauge 
invariant variables vanish, = $ = S' = in 7Z. 
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It is now straightforward to show that one can choose a gauge in which <f> = fa = 
7 = H = Si = A = 0. From the transformation laws (|A.5|I - i|A.ll|l we can choose 
C to make A = 0, choose A to make (j> = 0, and choose Bi to make fa = 0. The 
residual gauge freedom is then parameterized by functions A and Bi that are time- 
independent. Next, from Eq. (| A. 1 3|> together with the definitions 12.45fl - (|2.47() it 
follows that 

= 9 = iff (A.14) 
= $ = -2-y (A.15) 

= 3i = ~ii. (A.16) 

Thus the only remaining non-zero pieces of the metric other than the TT piece are 7 
and Si, and these are both time- independent. Finally we can use the residual gauge 
freedom given by time- independent functions A and Bi to set to zero 7 and £j, by 
Eqs. I|A.7|) and (|A.10(1 . [For this purpose A and Bi will vanish on dV, by Eqs. (|A.1(> 
and HA.4JI . so tp and r\i vanish.] 
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